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ABSTRACT 

Collisions are the core agent of planet formation. In this work, we derive an analytic description 
of the dynamical outcome for any collision between gravity-dominated bodies. We conduct high- 
resolution simulations of collisions between planetesimals; the results are used to isolate the effects 
of different impact parameters on collision outcome. During growth from planetesimals to planets, 
collision outcomes span multiple regimes: cratering, merging, disruption, super-catastrophic disrup- 
tion, and hit-and-run events. We derive equations (scaling laws) to demarcate the transition between 
collision regimes and to describe the size and velocity distributions of the post-collision bodies. The 
scaling laws are used to calculate maps of collision outcomes as a function of mass ratio, impact angle, 
and impact velocity, and we discuss the implications of the probability of each collision regime during 
planet formation. 

Collision outcomes are described in terms of the impact conditions and the catastrophic disruption 
criteria, Q* RD - the specific energy required to disperse half the total colliding mass. All planet 
formation and collisional evolution studies have assumed that catastrophic disruption follows pure 
energy scaling; however, we find that catastrophic disruption follows nearly pure momentum scaling. 
As a result, Q* RD is strongly dependent on the impact velocity and projectile-to-target mass ratio in 
addition to the total mass and impact angle. To account for the impact angle, we derive the interacting 
mass fraction of the projectile; the outcome of a collision is dependent on the kinetic energy of the 
interacting mass rather than the kinetic energy of the total mass. We also introduce a new material 
parameter, c*, that defines the catastrophic disruption criteria between equal-mass bodies in units 
of the specific gravitational binding energy. For a diverse range of planetesimal compositions and 
internal structures, c* has a value of 5±2; whereas for strengthless planets, we find c* = 1.9±0.3. We 
refer to the catastrophic disruption criteria for equal-mass bodies as the principal disruption curve, 
which is used as the reference value in the calculation of Q* RD for any collision scenario. The analytic 
collision model presented in this work will significantly improve the physics of collisions in numerical 
simulations of planet formation and collisional evolution. 
Subject headings: planets and satellites: formation, methods: numerical 



1. INTRODUCTION 

Planet formation is common and the number and di- 
versity of planets found increases almo st daily (e.g. 



Borucki et al. 2011 Howard et al. 2011). As a result, 
planet formation theory is a rapidly evolving area of 
research. At present, observations principally provide 
snapshots of either early protoplanetary disks or stable 
planetary systems. Little direct information is available 
to connect these two stages of planet formation, there- 
fore, numerical simulations are used to infer the details 
of possible intermediate stages. However, the diversity 
of extrasolar planetary systems continues to surprise ob- 
servers and theorists alike. 

A complete model of planet formation has eluded 
the astrophysics community because of both incomplete 
physics in numerical simulations and computational con- 
straints. In order to make the problem of planet for- 
mation more tractable, the process is often divided into 
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separate stages, which are then tackled in isolation. This 
method has had some success. For example, A-body sim- 
ulations show that large (~ 100 km) planetesimals may 
grow into protoplanets of about a lunar mass on mil- 
lion year time scales (e.g., Kokubo & Ida 2002). Other 
simulations, focusing on later stages ot planet formation, 
created a variety of stable planetary systems f rom initial 
distributions of protopla net size bodies (e.g., Chambers 



2001 Agnor et al. 



stable planets has 



19991. Recently, the distribution ot 



jeen investigated in population syn 



thesis models (e.g. , Mordasi ni et al.|2009 Ida fc Lin 2010 
Schlaufman et al.llzUlul lAlibert et al.llzUllD. However 



given the complexity of planet formation, it is unsurpris- 
ing that the predictions from the first population synthe- 
sis models have been overturned by the rapidly growing 
catalog of exoplanets ( Howard et al.|20"T0 2011 ). Hence, 
the diversity of the extrasolar planets is still unexplained. 

At the heart of the standard core-accretion model of 
planet formation is the growth of planetestimals. The 
evolution of planetesimals is dominated by a series of in- 
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dividual collisions with other pla netesimals (e.g., Beauge 
& Aarseth 1990; L issauer| |l993). The outcome ol each 
collision depends on the specific impact conditions: tar- 
get size, projectile size, impact parameter, impact ve- 
locity, and some internal properties ol the target and 
projectile, such as composition and strength. In the 
past, direct global simulations ol planetesimal evolution 
have assumed very simplified collision models. In TV- 
body simulations, terrestrial planet embryos were shown 
to grow easily from an annulus of large planete simals if 



the only out come of collisions is merging (e.g., Kokubo 
|fc Ida 2002). However, the computational demands ol 



such numerical methods did not permit lor the tracking 
of the very large numbers of bodies necessary to be able 
to include direct calculations of the erosion of planetesi- 
mals. 

Statistical methods are required to describe the full 
popula tion of bodies from dust s ize to planets. For ex- 
ample, Kenyon & Bromley ( 2009 ) conducted simulations 
that included fragmentation but still relied upon a sim- 
ple collision model. Specifically, their simulations did not 
fully account lor the effects of the mass ratio, impact ve- 
locity, or impact angle on the collision outcome. In order 
to overcome these simplifications some previous studies 
have employed a multi-scale approach that includes di- 
rect simulations ol collision o utcomes within a top-level 



simulation ol planet growth ( Leinhardt fc Richardson 
~ 2009l|Genda et al ||2011a[ ). EmT 



2005 Leinhardt et al. 



ever, multi-scale calculations significantly increase the 
computational requirements. In addition, the numerical 
methods employed in previous studies were only valid 
lor a specific impact velocit y reg ime. In the case of|Lein 



hardt fc Richardson| ( |2005[ ) and |Leinhardt et al.| (2009), 
the collision model assumed subsonic collisions ana could 
not be extended pas t oligarchic growth. In the case ol 



Genda et al. (2011a), the technique assumed strength- 
less bodies and cannot be used in the early phases ol 
planetesimal growth. 

A general description ol collision outcomes that spans 
the growth from dust to planets is required to build a 
sell-consistent model for planet formation. In previous 
work, the description of collision outcomes drew upon a 
combination of laboratory experiments and limited nu- 
merical simulations of collis ions between two plan etary- 
scale bodies (see review by Holsapple et al. 2002). Col- 
lision outcomes themselves are quite diverse, and sev- 
eral distinct collision regimes are encountered during 
planet formation: cratering, merging/accretion, frag- 
mentation/erosion, and hit-and-run encounters. 

Individual collision regimes have been described in 
quite varying detail. In the laboratory, the erosive 
regimes (cratering an d disruption) have been studied 
most comprehensively (iHolsapplell 19931 iHolsapple et al. 



20021; however, even these regimes lack a complete de 



scription ol the dependence on all impact parameters 
(particularly mass ratio and impact angle) . Recently, nu- 
merical studies ol collisions between sell-gravitating bod- 
ies ol similar size have identified new types of collision 
outcom es including hit-and-run and mantle-stripping 



events (|Agnor fc~A sphaug 2004 al | Asphaug et al. 2006 



Marcus et, a ljp09j |ii0i0b| |Leinhardt, et al.||!Wl(J| |As-| 
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phaug||2010[ Kok ubo &: Genda||2010[ |Benz et al.||2007| 
Genda et al.|2011b* ). Up to this point our understanding 



the diversity ol collision outcomes in planet formation 
codes. In addition, the transitions between regimes are 
not clearly demarcated in the literature. 

In the work reported here, we present a complete 
description ol collision outcomes for gravity-dominated 
bodies. Using a combination of published hydrocode 
and new and published A-body gravity code simulation 
results, we derive analytic equations to demarcate the 
transitions between collision regimes and the size and 
velocity distribution of the post-collision bodies. We de- 
scribe how these scaling laws can be used to increase 
the accuracy of numerical simulations of collisional evo- 
lutio n without sacrificing efficien cy. In a companion pa- 
per ( Stewart fc Leinhardt||2011 ), we apply these scaling 
laws to the end stage ol terrestrial planet formation by 
analyzing the range of collision outcomes from recent N- 
body simulations. 

This paper is organized as follows: <f2] summarizes 
the numerical method for the new A-body simula- 
tions. Spjderives a general catastrophic disruption scaling 
law. Then, we develop general scaling laws for the size 
and velocity distribution of fragments in the disruption 
regime. Q defines the super-catastrophic and hit-and- 
run regimes. [|5] presents the transition boundaries be- 
tween collision outcome regimes from our numerical sim- 
ulations and our analytic model. Sj6] discusses the range 
of applicability of our results, areas needing future work, 
and the implications of the scaling laws on aspects of 
planet formation. The Appendix summarizes the imple- 
mentation of the scaling laws in numerical simulations 
of planet formation and collisional evolution. Table |A.1| 
presents the definitions ol variables and annotations used 
in this work. 

2. NUMERICAL METHOD 

In this section, we describe the numerical method used 
in the new impact simulations presented in this work. 
Simulations ol relatively slow subsonic impacts were con- 
ducted using a TV-bod y code with finite-sized spherical 
particles, PKDGRAV |Stadel||2001| , which has been ex- 



tween small bodies (e.g., 


Leinhardt et al. 


2000 Michel 


et al.||2001| |Leinhardt & 


Richardson 


2002 


|Durda et al. 


2004; Lcinhardt k Richardson 2005; 


Leinhardt & Stew- 


art||2009| jLeinhardt et al\ 


2010p. 



Both the target and projectile are assumed to be rub- 
ble piles: gravitational aggregat es with no bulk tensile 
strength ( Richardson et aLp002 ) . The rubble pile parti- 



cles are bound together purely by sell-gravity. The parti- 
cles themselves are indestructible and have a fixed mass 
and radius (for cases without merging). The equations 
of motion of the particles are governed by gravity and 
inelastic collisions. The amount of energy lost in each 
particle-particle collision is parameterized through the 
normal and tangential coefficients of restitution. The 
rubble piles are created by placing particles randomly in 
a spherical cloud and allowing the cloud to gravitation- 
ally collapse with highly inelastic particle collisions. Ran- 
domizing the internal structure of the rubble piles avoids 
spurious collision results du e to crystalline structure of 
hexagonal close packing ( see Leinhardt et al.|2000 Lein- 



hardt fc Richardson|2002 1 . The crystalline structure can 



of these new regimes has not been sufficient to implement 



cause large uncertainties in collision outcomes for super- 
catastrophic events. 



Collision Scaling Laws 
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All simulations had a target with radius of 10 km, mass 
of 4.2 x 10 15 kg, bulk density of 1000 kg m~ 3 , and escape 
velocity of 7.5 m s _1 . The current study includes four 
projectile-to-target mass ratios (7), four impact angles 
(#), and a range of impact velocities spanning merging 
to super-catastrophic disruption. These results for a sin- 
gle size target body are used to derive scaling laws for any 
size body in the gravity regime. The target and projec- 
tile are initially separated by the sum of their respective 
radii to ensure that the impact angle of the impact is 
unchanged from the initial trajectory. 

In order to resolve the size distribution after the colli- 
sions, each body needs a relatively high number of par- 
ticles (iV targ - 10 4 , N p ~ 250 - 10 4 depending on the 
mass ratio). However, large numbers of particles are also 
time consuming to integrate, especially in a rubble pile 
configuration where there is a high frequency of particle- 
particle collisions. Each simulation here uses high resolu- 
tion with inelastic particle collisions to resolve the initial 
impact. Once the velocity field is well established, the 
particles are allowed to merge with one another. Thus, 
our method has both accuracy and efficiency, resolving 
the size distribution to small fragments and completing 
the simulations as quickly as possible. 

We considered the possible influence of the time of the 
transition from inelastic bouncing to perfect merging. 
Figure [l] presents a test case using a head-on catastrophic 
impact Detween equal-sized objects (mass ratio 7=1). 
Each body has ~ 10 4 particles; in the inelastic bouncing 
phase, each particle has a normal coefficient of restitu- 
tion, e n = 0.5, and a tangential coefficient of restitution, 
et = 1, consistent with field observati ons and friction ex - 
periments on rocky materials (e.g., Chau et al. 2002 ) . 
At a certain time, the colliding particles are allowed to 
merge, producing one particle with the same mass and 
bulk density as the two original particles. If merging is 
turned on too early, the mass of the largest remnant is 
overestimated (green and cyan line s) due to a geometric 



1000 



effect known as runaway merging (Leinhardt & Stewart 
20091. The results of this numerical test show that the 



mass distribution is stable if merging is turned on after 
50 steps, where 1 step is one minute in simulation time 
in the frame of the particles. However, we choose to be 
conservative and merge after 250 steps of inelastic bounc- 
ing in all of the new simulations presented in this paper. 
All simulations were run for at least 0.2 years, at which 
point the size distribution had stabilized and clumps of 
rubble pile fragments were easily identifiable. 

Previous studies using PKDGRAV did not have the 
numerical resolution to determine an accurate size or 
velocity distr ibution of the collisio nal remnants (e.g., 
- 10 3 , 



2000). In this work, we 



Leinhardt et al. 
present more extensive simulations at an order of magni- 
tude higher resolution (AT targ ~ 10 4 ). Note that N = 10 4 
is high resolution for A^-body simulations of colliding 
rubble-pile bodies with bouncing particles. We con- 
ducted several resolution tests and find that the random 
error on the mass of the largest remnant is a few percent 
of the total system mass. Hence, the super-catastrophic 
impacts, where the largest remnant mass is a few per- 
cent of the total system mass, have the highest error. 
We achieve excellent reproducibility of the slope of the 
size and velocity distributions with the nominal resolu- 
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Fig. 1. — Cumulative size distribution of collisional debris after 
the catastrophic impact between two 20-km diameter bodies. Line 
colors represent different handoff times from inelastic bouncing to 
perfect merging for the outcome of collisions between pairs of PKD- 
GRAV particles. Each step corresponds to 1 minute in simulation 
time. The same initial impact is used for all distributions shown: 
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tion compared to the higher resolution tests. 

Note that A^-body simulations are inherently higher 
resolution compared to smoothed particle hydrodynam- 
ics (SPH) simulations. Our simulations resolve over a 
decade in fragment size, comparable to SP H simulations 
using an order of magnitude more particles ( |Durda et al. 
2004|. Fragments of radius 0.5 km are considered the 



smallest usable fragments in these simulations. 

In the following sections, we also include published re 
suits of s ubsonic and supersonic collisions from previ 
ous work ( Leinhardt k. Stewart||2009 " l |Agnor k. As phaug 
2004alb| pVTarcus et al.pOOOl |2(H0b| |Durda et al.pOOT 
20071 1 Jutzi et al.|2010| |Benz ; fc Asphairg|19991 |Benz 2000 
Stewart & Leinhardt 2009 ; Korycansky fc Asphaug|2009 



Nesvorny et al.|2006 Benz et al.|2007[ ) . Studies of super 



sonic collisions utilize shock physics codes, which include 
the effects of irreversible shock deformation. For com- 
putational efficiency, the shock code is generally used to 
calculate only the early stages of an impact event; after 
a few times the shock wave crossing time, the amplitude 
of the shock decays to the point where further deforma- 
tion is negligible. After the hydrocode step, the gravita- 
tional reaccumulation stage of disruptive ev ents has been 
calculated directly using PKDGRAV (e.g., Leinhardt fc 
|Stewart||2009l iDurda et al.||2004l IMichel et al.||2002|) or 



indirectly by iteratively s olving tor the mass bound to 
the largest fragment (e.g ., Benz & Asphaug 1999 Benz 
2000) |Marcus et al.||2009| 



3. RESULTS: THE DISRUPTION REGIME 

In our model, the boundaries between collision out- 
come regimes are defined using the catastrophic disrup- 
tion criteria, the specific energy required to gravitation- 
ally disperse half the total mass, because it provides a 
convenient means of calculating the mass of the largest 
remnant. Our definition of the disruption regime refers 
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to collisions in which the energy of the event results in 
mass loss (fragmentation) between about 10% and 90% 
of the total mass. More quantitatively, the disruption 
regime is defined as collision that result in the largest 
remnant having a linear dependence on the specific im- 
pact energy. The rationale for this definition will become 
apparent in j |3.1.1| 

This section focuses on deriving the dynamical out- 
come (the mass and velocity distribution of post-collision 
fragments) in the disruption regime. Other collision 
regimes are discussed in Sj4j Before discussing the results 
of our new numerical simulations, we briefly review the 
catastrophic disruption criteria, as it is a fundamental 
part of our story. 

In the literature on planetary collisions, Q tradition- 
ally denotes the specific energy of the impact (kinetic en- 
ergy of the projectile/target mass) and Q* indicates the 
catastrophic disruption criteria, where the largest rem- 
nant has half the target mass. Upon recognition that 
gravitational dispersal was important, Q* s and Q* D de- 
noted the criteria for shattering in the strength regime 
and dispersal in the gravity regime, respectively. All of 
the previous definitions for Q* assumed that the projec- 
tile mass, M p , was much smaller than the target mass, 
M targ ; however, in several phases of planet formation it is 
expected that M p ~ M targ . Therefore, in previous work, 
we developed a disruption criteria in the center of mass 
reference frame in orde r to study collisions between com- 
parably sized bodies ( Stewart fe Leinhardt| [2009). The 
subscript R was added in the modification of the specific 
energy definition to denote reduced mass. The center of 
mass specific impact energy is given by 



Q R = (0.5M p l/ 2 p + 0.5M targ l/ 2 targ )/M t , 



-0.5fiV?/M toU 



(1) 



where Mt. 



*tot = Jw p -i- J"targ, M i g the reduced mass 
AfpMtarg/Aftot, Vi is the impact velocity, and V v and 
Ut alg are the speed of the projectile and target with re- 
spect to the center of mass, respectively. At exactly the 
catastrophic disruption threshold, 



CRD 



--0.5fiV* 2 /M tot , 



(2) 



where we explicitly define V* to be the critical impact 
velocity required to disperse half of the total mass for a 
specific impact scenario (total mass and mass ratio). 

The catastrophic disruption criteria is a strong func- 
tion of size with two components: a strength regime 
where the critical specific energy decreases with increas- 
ing size and a gravity regime where the critical specific 
energy increases with increasing size. The transition be- 
tween regimes occurs between a few 100-m and few-km 
radius, depending on the strength o f the bodies (see Fig- 



ure 2, Stewart & Leinhardt 20091. A general formula 



for Q* Bn as a fun ction of size was derived by Housen & 



Holsapple ( 1990 ) using 7r-scaling theory, 



!RD 



(sypi) 3p(0+3)/(20+3) ^^v( 3_2 *v*^ 2-3 ^ + 

q g { Pl Gf m i^V*( 2 - 3 «(3) 

where the first term represents the strength regime and 
the second the gravity regime. Rci is the spherical radius 
of the combined projectile and target masses at a density 
of pi = 1000 kg m -3 . The variable Rci was introduced 



by Stewart & Leinhardt ( 2009 ) in order to fit and com- 
pare the disruption criteria for collisions with different 
projectile-to-target mass ratios and to account for bod- 
ies with different bulk densities (e.g., rock and ice). G 
is the gravitational constant; q s and q g are dimension- 
less coefficients with values near 1. S is a measure of 
the material strength in units of Pa s 3 '^ +3 ^, and the re- 
maining variables, 4> and /S, are dimensionless material 
constants. is a measure of the strain-rate dependence 
of th e material strength with values ra nging from 6 to 9 
(e.g., |Housen fe Holsapple 1990 1999). p, is a measure 
of how energy and momentum from the projectile are 
coupled to the target; p is constrained to fall between 
1/3 for pure momentum scaling and 2/3 for pure energy 
scaling (Holsa pple fc Schmidt|1987 ). Note that the form 
of equation [3] assumes that the projectile and target have 
the same density. 

In the strength regime, the largest post-collision rem- 
nant is a mechanically intact fragment. The catas- 
trophic disruption criteria decreases with increasing tar- 
get size because more flaws grow and coalesce during 
the longer loading duratio n in larger impact events (e.g., 
Housen fc Holsapple|1999 |. In the gravity regime, disrup- 
tion requires both fracturing and gravitation al dispersal 
( TMelosh fc Ryan[1997||Benz fc Asphaug|1999[ ); hence the 
disruption criteria increases with increasing target size. 
In this regime, the largest remnant is a gravitational ag- 
gregate composed of smaller intact fragments. In both 
regimes, the disruption criteria increases with impact ve 
locity because more of the impact kinetic energy is dissi 



pated by shock def ormation at higher velocities ( Housen 



& Holsapple 1990). This work focuses on the gravity 
regime; the strength regime will be the subject of future 
studies. 

Both equations [2] and[3] are satisfied by collisions at ex- 
actly the catastrophic disruption threshold. The general 
formula for catastrophic disruption given by equation [3] 
describes a family of curves that depend on size, impact 
velocity, and material parameters (q g and jl). Most pre- 
vious work fit the material parameters in equation [3] to 
planetary bodies of a particular composition under var- 
i ous assumpt ions (e.g., a fixed impact velocity). Next 
( ^3. 1.2 3.1.4), we present a general method to calculate 
the values tor the disruption energy and critical impact 
velocity for specific impact scenarios and materials. 

3.1. Derivation of a general catastrophic disruption law 
in the gravity regime 

3.1.1. The Universal Law 

In previous work, using simulations of head-on im- 
pacts, we determined the value of Q* RD for a particular 
pair of planetary bodies by fitting the mass of the largest 
post-collision remnant, M\ r , as a function of the specific 
impact energy, Q R . The simulations held the projectile- 
to-targct mass ratio fixed and varied the impact velocity. 
For a wide range of target masses, projectile-to-target 
mass ratios, and critical impact velocities, we found that 
the mass of the largest remnant is approximated by a 
single linear relation, 



Mfc/Mtot = -0MQr/Q*rd - 1) + 0.5, 



(4) 



where Q* RD was fitted to be the specific energy such that 



M k = 0.5M tot (Stewart & Lemhardt||2009| iLeinhardt & 
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Fig. 2. — Schematic of the collision geometry. The target is 
stationary and the projectile is moving from right to left with speed 
Vi. The impact angle, 8, is defined at the time of first contact as the 
angle between the line connecting the centers of the two bodies and 
the normal to the projectile velocity vector. The impact parameter 
is b = sin 9. 




Fig. 3. — Normalized mass of the largest post-collision remnant 
versus normalized impact energy for all collisions in the disruption 
regime. The impact energy is scaled by the empirical catastrophic 
disruption criteria Q'^ D (Table The solid lines are the uni- 
versal law for the mass of the largest remnant (equation Hi; see 
text for discussion of 1:1 oblique impacts. The symbol denotes the 
projectile-to-target mass ratio, and the color denotes the impact 
parameter. 



Stewart|[2009[) . We found that a single slope agreed well 
with results from both laboratory experiments and nu- 
merical simuhttiojis^J^^ dimensional anal- 
ysis by Housen & Holsapple ( 1990 1 supports the linear- 
ity of the largest remnant mass with impact energy near 
the catastrophic disruption threshold. Hence, we refer 
to equation [4] as the "universal law" for the mass of the 
largest remnant. 

However, the most likely collision between two plane- 
tary bodies is not a head-on collisio n; a 45° impact angle 
is most probable ( Shoemaker|1962| ) . The impact param- 
eter is given by b = shift, where (J is the angle between 
the centers of the bodies and the velocity vector at the 
time of contact (Figure [2]). The impact parameter has 



a significant effect on the collision outcome because the 
energy of the projectile may not completely intersect the 
target when the impact is oblique. For example, in the 
collision geometry shown in Figure [2j the top of the pro- 
jectile does not directly hit the target (above the dotted 
line). As a result, a portion of the projectile may shear 
off and only the kinetic energy of the interacting fraction 
of the projectile will be involved in disrupting the target. 
Thus, a higher specific impact energy is required to reach 
the catastrophic disruption threshold for an oblique im- 
pact. 

The new PKDGRAV simulations conducted for this 
study were used to develop a generalized catastrophic 
disruption law, as previous work did not independently 
vary critical parameters. Table [T] presents the subset of 
the simulations discu ssed in detail below; for a complete 
listing see Table A.3| in the Appendix. The simulations 
are grouped by impact scenario: fixed mass ratio and 
impact angle. The value for the catastrophic disruption 
criteria, Q%d, is found by fitting a line to the mass of 
the largest remnant as a function of increasing impact 
energy in each group. The prime notation in the catas- 
trophic disruption criteria indicates an impact condition 
that may be oblique (b > 0). 

With the new data, we first consider how impact angle 
influences the universal law for the mass of the largest 
remnant. Figure [3] presents the normalized mass of the 
largest remnant versus normalized specific impact en- 
ergy. Our previous simulati ons (all at b = 0) are shown 
on the same universal law in |Stewart fc Leinhardt ( 2009). 
Note that comparable mass collisions with b > need to 
be considered carefully (offset for emphasis in Figure [3]). 
Such collisions transition from merging to an in elastic 
bouncing regime (called hit-and-run, discussed in S 3.1.2 ) 
before reaching the disruption regime. As a result, the 
mass of the largest remnant has a discontinuity between 
M to t and M targ with increasing impact energy. So, in the 
case of equal-mass collision^] only the fragments with 
Mir < Mtarg are fit by a line of slope -0.5. 

Our new results demonstrate that the same universal 
law for the mass of the largest remnant found for head-on 
collisions can be generalized to any impact angle: 

M h /M tot = -QMQr/Qrd - 1) + 0-5. (5) 

In detail, the mass of the largest fragment for a spe- 
cific subset of simulations may deviate slightly from the 
universal law (Figure [3]). Note that the deviations vary 
between subsets, with some results systematically sloped 
more steeply and others sloped more shallowly. The de- 
viations in Mir /M to t from equation [5] are about 10% for 
near-normal impacts (6 = 0.00 and o = 0.35) and some- 
what larger and more varied for highly oblique impacts. 

Overall, the universal law provides an excellent repre- 
sentation for mass of the largest remnant for all disrup- 
tive collisions in the gravity regime. As a result, we have 
chosen to use the range of impact energies that satisfy the 
universal law for the mass of the largest remnant as the 
technical definition of the disruption regime. At higher 
specific impact energies, the linear universal law breaks 

1 A robust fit requires several points between O.lMtot and Mtarg. 
A similar procedure was applied t o fit the Q'^ n for the b = 0.5, 
7 = 1 and 7 = 0.5 simulations from |Marcus et al.| ( |2010b| that are 
shown in Figure [4] 
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down in a transition to t he super-catastrophic regime 
(Qr/Qrd — 1-^' see [4.1 1. At lower specific impact en- 
ergies, the outcomes are merging or cratering (§51). Us- 
ing our definition, the disruption regime encompasses less 
than a factor of two in specific impact energy. The out- 
comes in the disruption regime span partial accretion of 
the projectile onto the target to partial erosion of the 
target body. 

Note that the derived values for Q'r D are strong func- 
tions of both the mass ratio and the impact parameter 
(Table [T]). The catastrophic disruption energy rises with 
smaller projectiles and larger impact parameters. |Benz| 



& Asphaug (fl999) investigated the effect of impact pa- 
rameter on the disruption criteria; however, their study 
fixed the impact velocity and varied the mass ratio of the 
bodies. Hence, the individual roles of the impact param- 
eter and mass ratio cannot be discerned from their data. 
In the next two sections, the influence of each factor is 
isolated and quantified. 

3.1.2. Dependence of disruption on impact angle and 
derivation of the interacting mass 

In order to describe the dependence of catastrophic dis- 
ruption on impact angle, we introduce two geometrical 
collision groups (Figure [2t: non-grazing - most of the 
projectile interacts with the target, and grazing - less 
than h alf the projectil e interacts with the target. Fol- 
lowing Asphaug (2010), the critical impact parameter, 



b C r\t — 



R 



R 



(6) 



is reached when the center of the projectile (radius r) is 
tangent to the surface of the target (radius R). Grazing 
impacts are defined to occur when b > i>crit- 

When considering a non-grazing impact scenario with 
a particular b and 7, the collision outcome transitions 
smoothly from merging to disruption as the impact veloc- 
ity increases. For grazing impacts, however, the collision 
outcome transitions abruptly from merging (Mj r ~ M to t) 
to hit-and-run (M\ T ~ M targ ) and then (less abruptly) to 
disruption (see <|5]) . Thus, only collision energies that re- 
sult in M\ t < Altarg should be used in the derivation of 
Q'r D in the grazing regime, as done for the 7 = 1 results 
shown in Figure [3] 

During oblique impact events, a significant fraction of 
the projectile may not actually interact with the target, 
particularly for comparable mass bodies. For gravity- 
dominated bodies, the projectile is decapitated and a 
portion of the mass misses the target entirely. As a result, 
only a fraction of the projectile's total kinetic energy is 
deposited in the target, and the impact velocity must 
increase to reach the catastrophic disruption threshold. 

Using a simple geometric model, we derive the fraction 
of the projectile mass that is estimated to be involved in 
the collision. First, we define I as the projected length 
of the projectile overlapping the target. As shown in 
Figure [2j 

l + B = R + r, (7) 

where B = (R+r) sinO. Placing the origin at the bottom 
of the projectile on the center line and the positive z-axis 
pointing to the top of the page, the estimated projectile 
mass involved in the collision, rni n t e ract, is determined by 



integrating cylinders of height dz and radius a from to 
I along the z-axis, 



p I na dz, 
la 



(8) 



where p is the bulk density of the projectile. The radius 
of each cylinder can be defined in terms of the radius of 
the projectile and the height from the origin, 



{r-zf 



Then, 



^interact = /0(ttW - (tt/3)1 ). 

Dividing by the total mass of the projectile, M p , 

^interact 3W I 



4r 3 



(9) 
(10) 

(11) 



Thus, a is the mass fraction of the projectile estimated 
to be involved in the collision (see Table [l] for the values 
of a in our simulations). The entire projectile interacts 
with the target when R > b(r + R) + r; then I < R and 
a = 1. 

In order to account for the effect of impact angle on 
Q'rdi we include the kinetic energy of only the interact- 
ing mass. The appropriate reduced mass is then 



Mo 



<xMpM( 



targ 



aM v + Mi 



(12) 



targ 



Now consider the difference between a head-on impact by 
a projectile of mass M p at V* and a head-on impact by 
a projectile of mass aM p . At the same impact velocity, 
the impact energies between the two cases differ by the 
ratio of the reduced masses, 



Q'r 



p_ 

/'a 



Qr, 



(13) 



Next, in order to conserve the effective specific impact 
energy, the impact velocity must increase with increasing 
impact angle such that 



V 



p_ 

Ma 



V* 



(14) 



However, the disruption criteria itself depends on the 
magnitude of the impact velocity (equation p|. In other 
words, when the effective projectile mass changes, the 
change in the impact velocity required for disruption 
varies by more than the factor presented in equation |14| 
Combining these two effects leads to the relationship bc- 
tween the oblique and head-on disruption energy for a 
fixed mass ratio collision, 



:rd 



A, 

Ma 

JL 
/'«> 




2-3ft 



'RD 



2-3/2/2 



Q 



RD- 



(15) 



By definition, the critical impact velocity for an oblique 
impact must satisfy equation [2] 



V 



I2Qrd m < 



(16) 



Collision Scaling Laws 
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TABLE 1 

Summary of parameters and results from selected PKDGRAV simulations (for full list of simulations see Appendix) 



Mp 


6 


Vj 


Mi, 


M slT 


P 


Qr 




a 


M lT /M tot 


M slT /M tot 




- 


m/s 


Mtot 


Mtot 


- 


J/kg 


J/kg 


- 


Predicted 


Predicted 


1.00 


0.00 


24 


0.76 


0.004 


4.0 


7.2 x 10 1 






0.67 


0.01 


1.00* 


0.00 


30 


0.50 


0.01 


3.2 


1.1 x 10 2 


1.1 x 10 2 


1 


0.49 


0.01 


1.00 


0.00 


35 


0.12 


0.05 


2.9 


1.5 x 10 2 






0.31 


0.02 


1.00 


0.35 


17 


0.48 


0.46 


4.8 


3.6 x 10 1 






0.44t 


0.44t 


1.00 


0.35 


30 


0.23 


0.20 


3.1 


1.1 x 10 2 


5.4 x 10 1 


0.72 


0.33t 


0.33t 


1.00 


0.35 


45 


0.02 


0.01 


3.7 


2.6 x 10 2 






o.nt 


o.nt 


1.00 


0.70 


12 


0.50 


0.49 


- 


1.8 x 10 1 






0.50t 


0.50t 


1.00* 


0.70 


80 


0.28 


0.27 


1.1 


8.0 x 10 2 


4.7 x 10 2 


0.22 


0.36t 


0.36t 


1.00 


0.70 


150 


0.04 


0.01 


3.2 


2.8 x 10 3 






super-cat 




1.00 


0.90 


20 


0.50 


0.50 


- 


5.0 X 10 1 






0.50t 


0.50t 


1.00 


0.90 


400 


0.35 


0.34 


4.6 


2.0 x 10 4 


1.5 x 10 4 


0.03 


0.39t 


0.39t 


1.00 


0.90 


600 


0.25 


0.25 


3.2 


4.5 x 10 4 






0.25t 


0.25t 


0.25 


0.00 


30 


0.69 


0.01 


3.7 


7.2 x 10 1 






0.72 


0.01 


0.25 1 "* 


0.00 


40 


0.40 


0.02 


4.4 


1.3 x 10 2 


1.3 x 10 2 


1 


0.51 


0.01 


0.25 


0.00 


50 


0.09 


0.02 


3.4 


2.0 x 10 2 






0.23 


0.02 


0.25 


0.35 


30 


0.67 


0.01 


4.5 


7.2 x 10 1 






0.81 


0.005 


0.25* 


0.35 


40 


0.53 


0.01 


4.1 


1.3 x 10 2 


1.9 x 10 2 


0.93 


0.66 


0.01 


0.25 


0.35 


60 


0.25 


0.01 


3.8 


2.9 x 10 2 






0.24 


0.02 


0.25 


0.70 


50 


0.69 


0.01 


2.8 


2.0 x 10 2 






0.90 


0.003 


0.25 


0.70 


100 


0.52 


0.004 


3.3 


8.0 x 10 2 


9.9 x 10 2 


0.33 


0.60 


0.01 


0.25 


0.70 


150 


0.32 


0.01 


2.5 


1.8 x 10 3 






0.09 


0.02 


0.25 


0.90 


120 


0.77 


0.14 


3.17 


1.2 x 10 3 






0.94** 


0.002** 


0.25 


0.90 


350 


0.47 


0.003 


4.40 


9.9 x 10 3 


9.3 x 10 3 


0.05 


0.47 


0.01 


0.25 


0.90 


450 


0.31 


0.01 


3.28 


1.6 x 10 4 






0.13 


0.02 


0.10 


0.00 


40 


0.79 


0.001 


4.9 


6.7 x 10 1 






0.79 


0.005 


0.10 


0.00 


65 


0.41 


0.01 


3.7 


1.8 x 10 2 


1.6 x 10 2 


1 


0.45 


0.01 


0.10 


0.00 


80 


0.14 


0.03 


3.4 


2.7 x 10 2 






0.17 


0.02 


0.10 


0.35 


40 


0.79 


0.002 


3.7 


6.7 x 10 1 






0.88 


0.003 


0.10 


0.35 


80 


0.47 


0.01 


4.5 


2.7 x 10 2 


2.7 x 10 2 


ltt 


0.51 


0.01 


o.iot* 


0.35 


100 


0.33 


0.01 


3.6 


4.2 x 10 2 






0.23 


0.02 


0.10* 


0.70 


100 


0.77 


0.002 


3.6 


4.2 x 10 2 






0.90 


0.003 


0.10 


0.70 


200 


0.52 


0.004 


3.8 


1.7 x 10 3 


2.0 x 10 3 


0.46 


0.59 


0.01 


0.10 


0.70 


300 


0.21 


0.01 


3.3 


3.7 x 10 3 






0.07 


0.02 


0.10 


0.90 


400 


0.70 


0.001 


5.2 


6.7 x 10 3 






0.89 


0.003 


0.10 


0.90 


700 


0.53 


0.002 


4.8 


2.0 x 10 4 


2.9 x 10 4 


0.07 


0.65 


0.01 


0.10 


0.90 


900 


0.57 


0.01 


4.3 


3.4 x 10 4 






0.42 


0.02 


0.025 


0.00 


100 


0.77 


0.001 




1.2 x 10 2 






0.91 


0.002 


0.025 


0.00 


140 


0.55 


0.01 


4.45 


2.3 x 10 2 


6.4 x 10 2 


1 


0.82 


0.01 


0.025 


0.00 


160 


0.51 


0.01 


4.10 


3.1 x 10 2 






0.76 


0.01 


0.025 


0.35 


160 


0.60 


0.01 


3.97 


3.1 x 10 2 






0.79 


0.01 


0.025 


0.35 


200 


0.45 


0.01 


4.78 


4.8 x 10 2 


7.2 x 10 2 


ltt 


0.67 


0.01 


0.025 


0.35 


300 


0.07 


0.05 


3.25 


1.1 x 10 3 






0.25 


0.02 


0.025 1 "* 


0.70 


300 


0.65 


0.002 


5.08 


1.1 x 10 3 






0.73 


0.01 


0.025 


0.70 


400 


0.47 


0.01 


3.59 


1.9 x 10 3 


2.0 x 10 3 


0.74 


0.52 


0.01 


0.025 


0.70 


500 


0.26 


0.01 


3.23 


3.0 x 10 3 






0.26 


0.02 


0.025* 


0.90 


800 


0.74 


0.001 




7.7 x 10 3 






0.65 


0.01 


0.025t* 


0.90 


900 


0.66 


0.002 


4.06 


9.7 x 10 3 


1.1 x 10 4 


0.12 


0.56 


0.01 


0.025 


0.90 


1000 


0.36 


0.003 


3.77 


1.2 x 10 4 






0.46 


0.01 



^ p - mass of projectile normalized by mass of target; b - impact parameter; Vi — projectile impact velocity; j^ lr _ mass of largest 

remnant normalized by total mass; M s \ r - mass of the second largest remnant; ft - slope of cummulative size distribution; Qr - center of 
mass specific energy; Qr D - empirical critical center of mass specific energy for catastrophic disruption and gravitational dispersal derived 
from the simulations. In all cases, the target contained ~ 1 X 10 4 particles, M targ = 4.2 X 10 15 kg, iJt ar g = 10 4 m; * indicates models 
shown in blue in Figure [H] t indicates N\ T = 2 and N s \ r = 4; ** erosive hit-and-run regime, the disruption regime model does not apply; 
tt indicates an a for which b > but ( < _R thus o = lj — not enough material to fit a power law; t* indicates models shown in Figurer?] 
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The correction for changing the mass ratio is derived in 
the next section. 

Our model for the effect of impact angle is used to 
derive equivalent head-on Q* RD values from our new and 
previously published catastrophic disruption data. Using 

the values for Q RD and V * fitted to the oblique simu- 
lation results, the equivalent head-on impact disruption 
criteria are 



Q 



RD 



v* 



I RD 



(3/1/2-2) 



2 Q*RD M t' 
P 



(17) 



(18) 



We considered the catastrophic disruption of a wide 
variety of planetary bodies from the studies summarized 
in Table |A~2j First, we fit the general expression for Q* RD 
(equation 1 3; n) to the (equivalent) head-on disruption data 
to derive the values of q g and p that best describe the 
entire data set, from planetesimals to planets. The same 
value for the material parameter /2 is used in the angle 
correction and the fit to equations] A small number of 
data points were ex cluded from the global fit, which are 
discussed in { 6.2.1 The best fit values for q g and p were 
found by minimizing the absolute value of the log of the 
fractional error, S = \\og(Q RDsi JQ RDmodcl )\. 

In some cases, the impact angle correction is signifi- 
cant (e.g., the impact scenarios with small values of a 
given in Table [I]) . With the exception of the con stant- 
velocity resu lts from Benz & Asphaug (|1999[) and Jutzi 
et al. ( 20T0J ) and the mixed velocity data from ( jfienz 
2U(Jl)p , the disruption data were derived from simulations 
conducted with a constant mass ratio and the critical im- 
pact velocity for catastrophic disruption, V * , was found 
by fitting to the universal law. For the simulations de- 
scribed in Table [l] the model correction for impact angle 
usually yields an impact energy within a factor of 2 of 
the simulation results for head-on collisions (e.g., within 
the linear regime for the mass of the largest remnant). 
We restricted our fits to cases where a > 0.5 to reduce 
any error contribution from a poor model correction for 
highly oblique impacts. 

The compiled data and best fit model Q* RD are pre- 
sented in Figure and B. The combined data are well 
fit by q g = 1.0 and p = 0.35 (5 = 0.14). Note the good 
match in the values for V* (colors) from the simulations 
with the lines of constant V*. Similarly good fits are 
found for 0.33 < p < 0.36 and 0.8 < q g < 1.2 with 
0.14 < 5 < 0.15. Amazingly, the compilation of catas- 
trophic disruption data is well fit by equation [3] for single 
values of q g and p for a wide variety of target composi- 
tions and over almost 5 orders of magnitude in size and 
9 orders of magnitude in impact energy. The critical im- 
pact velocities span 1 m s _1 to several 10's km s _1 . The 
best fit value for p falls near pure momentum scaling 
(/i=l/3). 

Upon closer examination, we found that the global fit 
with equation [3] systematically predicts a low disrup- 

2 In the fitting procedure, the strength term is neglected in equa- 
tion [3] For the lines plotted in Figure [4] the strength regime pa- 



rameters ar e fixed at <j> = 7, S = 2.4 Pa s 0,3 
the work in |Stewart fc Lei nhardt ( 2 009 \ . 



and q s 



1 based on 



tion energy for small bodies (Rci < 1000 km) and a 
high disruption energy for planet-sized bodies. Next, we 
consider separately the data for small and large bod- 
ies. A better fit is found for the small body data in 
Figure HA with 0.35 < p < 0.37 and 1.4 < q g < 1.65 
with 0.11 < S < 0.12. The small body data includes 
hydrodynamic to strong bodies and different composi- 
tions. The planet data in Figure [4j\ are best fit with 
0.35 < p < 0.375 and 0.85 < q g < 1.0 with the very 
small error of 0.038 < 5 < 0.041. The planet size data in- 
cludes three different target compositions. The data for 
collisions between small strong bodies ha ve the largest 
dispersion; these data will be discussed in §6.2. 1| 

3.1.3. Dependence of disruption on mass ratio 

By fitting such a large collection of data, it is clear that 
equation [3] describes a self-consistent family of possible 
Qrd values. For a specific impact scenario, the correct 
value for V* at each Rci is ambiguous because V* de- 
pends on both a material property and the mass ratio. In 
studies that hold Vt constant and vary the mass ratio, the 
derived value for the critical impact energy only applies 
for the corresponding critical mass ratio. As noted in pre- 
vious work, the critical impact vel ocity falls dramatically 
as the mass rati o approaches 1:1 ( Benz pOOO Stewart & 
Leinhardt 20091. As a result, collisions between equal- 
mass bodies require the smallest impact velocity to reach 
the catastrophic disruption threshold. 

Because a mass ratio of 1:1 defines the lowest disrup- 
tion energy for a fixed total mass, we derive the disrup- 
tion criteria for different mass ratios with respect to the 
equal- mass disruption criteria, Q* RD 7=1 . We begin with 
the equality between the impact energy and gravity term 
in the disruption energy (equation [3]) , 



Qr — Qrdi 



2M t , 



q g ( Pl G) 3fi/2 Rf^V-W . 



Note that 



and 



H = M p M targ /(M p + M targ ), 
7 



7 + r 



-Mi 



targ? 



M tot = (7 + l)M targ 
Then, substituting for // and M to t, 
( 7 /(7 + l))M targ U* 2 



(19) 

(20) 
(21) 



2( 7 + l)M t 



3/^/2 R 3/x T/ *(2-3/2) 



-q g ( Pl G)WR%[V 



arg 



v* 



2(7+ 1) 2 

7 

1 (7 + I) 2 

4 7 



q g ( Pl G)^ /2 Rf 1 



1 1/(3/*) 



1/(3/2) 



V* i 



(22) 



Then, for the same total mass, the relationship between 
the equal-mass disruption energy and any other mass 
ratio is determined by the difference in the critical impact 
velocities, 



Qrd — Qrd,j=i 



V* 



(2-3/2) 



Collision Scaling Laws 
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Fig. 4. — A compilation of gravity-regime catastrophic disruption simulation results. Symbols denote different target materials (Table |A.2|l , 
and color denotes the critical impact velocity, V* . Filled and line symbols are head-on impacts; open sy mbo ls a re o blique impacts" A 
and B: Data corrected for impact angle to equivalent head-on impact using the interacting mass (equat ions 1 1 7| and 1 1 8 | l . Constant-velocity 
Q*RD curves (equation pi are best fit to all the data with p, = 0.35 and q g = 1. The fit between the data and model curves is very good 
over almost 5 orders oimagnitude in size and 9 orders of magnitude in impact energy. Contours for V* = .005, .02, .1, .3, 1.5 and 5 
km s _1 (A k, C) and V = 15, 20, 30, 40, 60, 80 km s" 1 (B & D). C and D: Data converted to an equivalent equal -mass (1:1) disruption 
criteria using equations |23| and |22| The equal-mass data fall on lines proportional to Rq 1 - Fits to the equal-mass data are called "principal 
disruption curves" that arc denned by c* (black lines, equation [28[l; c* represents the value for the equal-mass Q% D in units of the specific 
gravitational binding energy. Best fit values are c* = 5 and p, =HX37 for small bodies and c* = 1.9 and p, = 0.36 for hydrodynamic planets. 
Inset: Full Q* RD curves (0.1, 1, 10, 100 km s — 1 ) showing transition from strength to gravity regimes. 

2/(3j2)-l 



= Q 



RD,~/=1 



1 (7 + I) 2 

4 7 



(23) 



tion. 



3.1.4. 



The equations for Q* RD 7=1 and V*—-y are given in the 
next section. 

In the compilation of catastrophic disruption data 
shown in Figure |4p and D, all the 7 < 1 data have 
been converted to an equivalent equal-mass impact dis- 
ruption energy and the colors denote V* =1 . For example, 
the critical disruption energy from head-on PKDGRAV 
simulations with 7 = 0.03 are a factor of three above 
the disruption energy for 7 = 1 in Figur e HK (parallel 
sets of + from Stewart fc Leinhardt|2009| . The data lie 
on the same line atter the correction in Figure [4p. The 
correction also brings together data from studies using 
different numerical methods and vastly different mate- 
rial properties. For example, the high-velocity Q* RD for 
strong and weak basalt targets (Y < ®M) fall on the same 
line as the PKDGRAV rubble piles after the conversion 
to an equivalent equal-mass i mpact. Similarly, studies 
of the disruption of Mercury (Benz et al. |2007 ) follow 
the sam e curve as disruption of earth-mass water/rock 
planets ( Marcus et al.pOlOb ). The general form for the 
equal-mass disruption criteria is derived in the next sec- 



The principal disruption curve 

In the previous two sections, we calculated the disrup- 
tion criteria for head-on equal-mass collisions by adjust- 
ing the critical disruption energy to account for different 
impact angles and mass ratios. The head-on equal-mass 
data points, derived from the compilation of numerical 
simulations, fall along a single curve that we name the 
"principal disruption curve" (black lines in Figure Wp 
and D). 

On the principal disruption curve, the critical impact 
velocity for equal-mass head-on impacts, V* =1 , satisfies 
both equation [l] and the gravity regime term in equation 



m 



^7=1^7=1 



'.R-D,7=l 



2M 



tot 



(24) 



Then, substituting \i 1=1 = M targ /2 = M tot /4, 



V* 



8g g ( P1 G) 



3/1/2 



7? 3ai 
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1/2 



(25) 



Thus, along a curve with a fixed projectile-to-target mass 
ratio, the critical impact velocity has a linear dependence 
on Rci ■ The linear dependence of V* on Rci for a fixed 
m ass ratio was confirmed by the numerical simulations 
in jStewart fe Leinhardt] ( |2009| ) (+ in Figure [§. 



_'hen, consider the dependence of the catastrophic dis- 
ruption criteria on size (equation [3]) and replace the ve- 
locity term with size, 



ex. n, ci n ci 
ocR 2 



m 

1 

(2-3/i) 



CI- 



(26) 

Thus, the catastrophic disruption criteria scales as radius 
squared along any curve with a fixed projectile-to-target 
mass ratio. 

Next, note the proximity of the gravity-regime equal- 
mass disruption energy to the specific gravitational bind- 
ing energy, 

u = ^t- (27) 

shown as the grey line in Figure |4j We define a dimen- 
sionless material parameter, c*, that represents the offset 
between the gravitational binding energy and the equal- 
mass disruption criteria. Then, the principal disruption 
curve is given by 



!RD,j=l 



4 

'-npiGR, 
5 



■ 2 

ci- 



(28) 



The parameter c* is a measure of the dissipation of en- 
ergy within the target. 

The coefficient q g is found by substituting equation [2] 
for Q* RD 7=1 in equation 28 and then equation [25] for 



V* 



2M tot 



'-tt Pi GR c1 , 



A 

/327TC* 



3m/2 



Finally, substituting q g into equation 25 gives 



V* , = 



/32 7 



V 5 



1/2 



(p 1 G) 1 / 2 i? C i 



(29) 



(30) 



Hence, the critical velocity along the disruption curve for 
equal-mass impacts is solely a function of Rci and c*. 

The principal disruption curve (equation 28 ) is a sim- 
ple, yet powerful way to compare the impact energies 
required to disrupt targets composed of different mate- 
rials. Each material is defined by a single parameter c*. 
In Figure |4p and D, the best fit values are c* = 5 ± 2 
and \x — 0.37 ± 0.01 for small bodies with a wide va- 
riety of material characteristics and c* = 1.9 ± 0.3 and 
p, = 0.36 ± 0.01 for the hydrodynamic planet-size bod- 
ies. These simulations span pure hydrodynamic targets 
(no strength), rubble piles, ice, and strong rock targets. 
Hence, for all the types of bodies encountered during 



planet formation, c* is limited to a small range of values. 
Note that the difference in c* between the small and large 
bodies is not simply because of the differentiate d struc- 
ture of the large bodies; two pure rock cases (►, Marcus 
et al. 2009) fall on the same Q*rj) curve. Kather, 
the large bodies were all studied using a pure hydro- 
dynamic model, whereas the small bodies were studied 
using techniques that incorporated material strength in 
various ways. A transition from a higher value for c* for 
small bodies to a lower value for planet-sized bodies is 
appropriate for planet formation studies, as discussed in 

m 

Now it is clear that most of the differences in the catas- 
trophic disruption threshold found in previous work are 
the result of differences in impact velocity and mass ratio 
(few studies varied impact parameter). 

Here, we have derived a general formulation for the 
catastrophic disruption criteria that accounts for mate- 
rial properties, impact velocity, mass ratio, and impact 
angle. The forward calculation of Q' RD f° r a specific im- 
pact scenario between bodies with material parameters 
c* and jl is de scribed in the Appendix an d in the com- 
panion paper ( Stewart k. Leinhardt||20ilj ). 



3.2. Fragment size distribution 

In the disruption regime, our new simulations resolve 
the size distribution of fragments over a decade in size 
(Figure [5]). In general, the post-collision fragments 
smaller than the largest remnant form a smooth tail 
that can be fit well by a single power law. The second- 
largest remnant forms the base of this tail. For most 
collisions there is a significant separation in size between 
the largest and second largest remnant. However, if the 
collision is very energetic, the largest remnant joins the 
power-law distribution (e.g., in 7 = 1,6 = 0.35). In 
addition, for the hit-and-run impacts with 7=1, the 
two largest remnants are comparable in size. Only the 
most energetic scenarios with 7 = 1 fall in the disrup- 
tion regime (e.g., black lines in b — 0.35 and 0.7). In all 
disruption regime collisions, the slope of the cumulative 
power-law tail, —/3 (see Tabic HI), is effectively indepen- 



dent of the impact conditio ns (0, Vj, 7) 



Using the method from |Wyatt fc Deht ( |2002[ ) and 
Paardekooper, Leinhardt, and Thebault (in prepara- 
tion), the mass of the second largest remnant, M s i T , is 
fully constrained by knowledge of the mass/size of the 
largest remnant and the power-law slope for the size dis- 
tribution of the smaller fragments. Let us consider a 
differential size distribution 



i(D)dD = CD- {fi+1) dD, 



(31) 



where n(D) is the number of objects with a radius be- 
tween D and D + dD, —((3 + 1) is the slope of the dif- 
ferential size distribution, and C is the proportionality 
constant. Integrating equation |31| the number of bodies 
between D\ T and D S \ T is 



N(D h ,D sh .) = -^(Dj-D-/ 



(32) 



Therefore, the number of objects larger than the second 
largest remnant (£> s ir) is N s \ r = N(D S \ T , 00). Assuming 
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Fig. 5. — Cumulative size distribution versus fragment diameter. For each mass ratio 7 and impact parameter 6, size distributions are 
shown for three different impact energies. Table ^provides the details for these simulations. The colors are an aid for the eye: magenta is 
the lowest energy impact in each block of three in Table [T] black is the highest energy, and cyan is in-between. In five panels, the fragment 
size distribution scaling law (blue line and triangles) is compared to the data. The impact parameters used for the model comparison are 
indicated in Table ^ by a * in the first column. 



that /3 > 0, 



N sh 
C 



(33) 



For spherical bodies with bulk density p, the mass of 
material between D\ T and D S \ T is 



M(D h ,D sh .) 



r pc 3-/3 • 



(34) 



In order to enforce a negative slope of the remnants, 
j3 must be less than 3. Mass is conserved in the im- 
pact; thus, the mass in the remnant tail must equal the 
total mass minus the mass in the largest remnant (s), 
M(0, .D s i r ) = M tot — N\ T M\ T , where N\ v is the number 



of objects with mass equal to the largest remnant (here, 
we allow for multip le l argest remnants) . Substituting for 
I? s i r from equation |33[, C is given by 



C 



N sh f3 



(3-^)(M to t-iVlrMi r ) 

(4/3)irpN sll /3 



(35) 



Substituting this expression for C into equation 33 and 
assuming that all of the objects are spherical, the size 
and mass of the second largest remnant is expressed in 
terms of the total mass by 



D S \r 
D tn t 



(3-0X1 



iy lr M > 



N sh (3 



(36) 
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where Dtot = 2((3M to t)/(47rp)) 1 / 3 . In the simulations 
presented here, the calculated diameter of the frag- 
ments is not informative because most PKDGRAV par- 
ticles merge with other particles in gravitationally bound 
clumps; in these cases, the bulk density is assumed for 
the size of the merged particle. The mass of the remnants 
is accurate, however. Rewriting equation [36] in terms of 
mass, 



M 



sir 



(3-/3)(l-iV lr ^) 



N sh (3 



(37) 



where M\ r is given by the universal law and the catas- 
trophic disruption criteria Q^qfequation 5]). 

In the last column of Table [T] thepredicted mass of 
the second largest remnant (equation 37 ) is compared to 
the numerical simulations using the empirically fit Q'r D , 
f3 = 2.85, Nit = 1, and N s \ r = 2. Since the analytic 
method presented here assumes an infinite size distri- 
bution in the fragment tail, we selected the value of /3 
to optimise the fit to the value of M s \ r in the simu- 
lations. This simple method of predicting M s i r works 
well for all impact conditions^] To illustrate the model, 
the predicted size distribution (blue line and triangles) 
is compared to selected numerical simulation results in 
Figure [5] In order to predict the size distribution of frag- 
ments in the hit-and-run regime impacts between com- 
parable mass bodies (7 = 1 and b > 0), the model needs 
to be modified slightly. In this special case, we suggest 
adopting Ai r = 2 and N s i T = 4 because the target and 
projectile each have a nearly identical size distribution 
of fragments. In this example, the model is calculated 
for the same impact conditions as the cyan data s et w ith 



4.2 



for 



7 = 1, b = 0.7, and Vi = 80 m s" 1 (see section 
more detailed discussion of the hit-and-run regime)? 

The fragment size distributions calculated using our 
subsonic TV-body simulations are consistent with shock 
code calculations investigating as teroid family formation 



via catastrophic i mpact events (Nesvorny et al. 2006 
All of the asteroid family-forming 



Jutzi et al. 



2010) 



simulations used a hybridized numerical technique, com- 
bining an SPH code with PKDGRAV in order to model 
the propagation of the initial shock wave and the subse- 
quent gravitational reaccumulation of the collision rem- 
nants. The asteroid family-forming collisions have sig- 
nificantly different impact parameters compared to our 
simulations: Vi was orders of magnitude larger, 7 was an 
order of magnitude smaller than our smallest 7, and tar- 
gets were larger (10's km in diameter). These differences 
notwithstanding, we find that the range in the values 
of j3 for the tail of the size distribution is very similar 
to our A-body results (note that some published values 
for /? include the largest remnant in the fit, whereas we 
do not). Qualitatively, we also find a general t rend in 
curvature of the size distribution consistent with IDurdal 



3 Because the analytic model for the fragment size distribution! 
assumes laniniimTeraiigetrfsiz^ 



be less than 3, which is slightly smaller than the slope of power 
laws that are best fit to the data ('fable [Til. The fragment size 
distribution may be modeled under different assumpti ons : 
choosing a minimum diame ter in the integral of e quation^ 
would represen t the smallest co nstituen t particles or grain 
We have chos en not the impose any ass umptions about material 
properties in the model presented here, but there may b e situations 
were more is known about the colliding bodies and the model for| 
determining A! s i r and fi may be modified. 
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Fig. 6. — Velocity of largest remnant with respect to the initial 
center of mass target velocity versus the mass of largest remnant 
normalized by the mass of the target. Impact angle is indicated by 
color; mass ratio is indicated by symbol. 



et al. (20071, with slightly convex size di strib utions for 
super-catastrophic impact events (section 4.1). 



3.3. Fragment velocity distribution 

Next, we consider the velocity of the collision rem- 
nants. The results are easier to interpret by separating 
the velocities of the largest remnant from the rest of the 
collision remnants. 

We first consider the speed of the largest remnant with 
respect to the center of mass of the collision (Figure [6]). 
For erosive events (M\ T < M targ ), there is almost no 
change in the amplitude of the target velocity for impacts 
with b = 0.9. Even at b = 0.7, the velocity reduction is 
minimal for all fractions of mass lost. Because b = 0.7 is 
the center of the probability distribution of impact an- 
gles, fully half of all erosive impacts have < 10% change 
in the target velocity amplitude. After head-on colli- 
sions (6 = 0), the largest remnant moves with the center 
of mass velocity. Note there is significant scatter in the 
data from 7 = 0.025, which is due to the fact that there 
was a small number of particles delivering the impact 
energy to a localized region of the target; thus, the orga- 
nization of the surface features on both objects become 
important. For disruptive impacts at b = 0.35, there is 
partial velocity reduction of the largest remnant. From 
these data, we cannot define a unique function for the de- 
pendence of V\ r on &, and we suggest that a quasi-linear 
relationship for < b < 0.7 is a reasonable approxima- 
tion. We stress that the specific dependence of Vi r on b 
in the disruptive regime is likely to be sensitive to inter- 
nal structure and composition, so extrapolation of these 
results beyond weak, constant density objects should be 
done with caution. 

In complete merging events, of course, the post-impact 
velocity is zero with respect to the center of mass. The 
b = 0.35 data with M\ r > Mt arg steadily approach the 
center of mass velocity with more mass accreted. The 
b = 0.7 and 0.9 data points plotted near Mi r /M targ = 1 
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Fig. 7. — Fragment mass— velocity histograms for simulations in Figure [5] and Table [T] The fragment velocities are relative to the largest 
remnant in units of the escape velocity from the combined mass of the target and projectile, V csc = (2GMtot/ Rci) 1 ^ 2 ■ The color coding 
is the same as in Figure [5] The scaling law predictions are shown in blue. 



are prim arily hit-and-run events, which will be discussed 
in j p~2l 

The smaller remnants of disruptive collisions have a 
more complex behavior. Figure [7j presents mass his- 
tograms of fragments versus velocity with respect to the 
largest remnant from the simulations summarized in Ta- 
ble [l] The slowest simulations are not plotted for the 
7 = 0.25 and 1.0 grazing impacts because there are only 
a small number of fragments. A significant number of 
the fragments consist of 10 PKDGRAV particles or less; 
in Figure [7J mass associated 10 or less particles is shown 
as dotted histograms. The dotted histograms overlay the 
total mass histograms for all but the lowest velocity bins; 
thus, within most of the velocity bins, the simulations 
do not have the resolution to robustly constrain the size- 
frequncy distribution of the mass in the bin. The smallest 
(poorly resolved) fragments are found in all velocity bins, 



while the largest fragments tend to move slowly with re- 
spect to the largest remnant. For example, the second 
largest remnant falls in one of the lowest velocity bins, 
but that bin is also occupied by smaller fragments. 

Hence, to describe the velocity field after a collision, 
we fit the velocity-binned mass of the collision remnants. 
The binned mass versus velocity is a fairly well defined 
exponential function for most of the simulations. In 
general, the lowest velocity bin in Figure [jj is of order 
0.1M tot . Using a least-squares fit of the subset of simu- 
lations in Table [TJ we find the mass fraction in the lowest 
velocity bin is proportional to the largest remnant mass: 



A = -0.3M lT /M u 



0.3. 



(38) 



To determine the slope, S, of the binned mass versus ve- 
locity exponential function, we integrate the differential 
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Fig. 8. — The mass of the largest remnant in the catastrophic 
and super-catastrophic disruption regimes. The solid line shows 
the combined universal law (equation™ an d re commended power 
law relation for M\ r /Mt a t < 0.1 (equation |44[ l. The symbols are 
new gravity regime simulations and the dotted and dashed lines 
represent the range of super-catastrophic disruption data in labo- 
ratory experiments in the strength regime. The shape and color of 
the symbols are the same as in Figure [3] 

mass function, 



a dm 
Av — = 10 
dv 



A-Sv 

5 

dm 10 A ' Sv 



dv 
M rem _ 
M tot " 

S-- 



Av 







A 



dv, 



v 



ln(10) Av (M rem /M tot ) ' 



(39) 
(40) 
(41) 
(42) 
(43) 



39] is 

icated 



where m = M/M to t, v = V/V csc , Av is the bin width, 
and the total mass in the histogram is the total mass in 
the remaining remnants, M rem = Af to t — M\ T . 
The fragment velocity scaling law (equation 

shown in blue in Figure [7] for selected cases inc 

by t* in Table [ll The velocity distributions of the rem 
nants agree qualitatively with those found in hyperve- 
locity simulations of asteroid family forming events, al- 
though previous workers have not fit any fu nction to the 
veloc i ty distribution of the fragments (e.g., Michel et al. 
2002| |Nesvorny et al.|2006|. 



4. OTHER COLLISION REGIMES 
4.1. Super- catastrophic regime 
In both laboratory experiments in the streng th regime 



(e.g., Kato et al. 1995 Matsui et al. 19821 and the 



few high resolut ion disruption simul ations in the grav- 



ity regime (e.g., Korycansky & Asphaug 



20091, the re- 



lationship between the mass of the largest remnant and 
the specific impact energy Qr shows a marked change in 



slope at around M\ r /M to t ~ 0.1. We define the super- 
catastrophic regime when Mi r /M tot < 0.1 (e^g-, when 
Qr/Q'rd > 1-8 by the universal law, equation pM) . In the 
super-catastrophic regime, the mass of the largest rem- 
nant follows a power law with Qr rather than the linear 
universal law. 

The slope of the power law for the largest remnant 
mass vs. impact energy shows some scatter in laboratory 
data, primarily in the range of -1.2 to -1.5. In Figure [HJ 
our few simulations of super-catastrophic collisions (sym- 
bols) are compared to the range of outcomes from labo- 
ratory experiments (dotted and dashed lines). Based on 
the simulations in the gravity regime and laboratory ex- 
periments in the strength regime, we recommend a power 
law in the super-catastrophic regime, 

M h /M tot = ^(Qr/Q'Zd)*, (44) 

where r\ *~ —1.5 and the coefficient is chosen for conti- 
nuity with the universal law (equation [5| . The slope of 
the power law, about -1.5, is consistent with our grav- 
ity regime simulations and a wide ran ge of laboratory 
experi ments summarized in Figure 1 in |Holsapple et al. 

pool 

In Figure [8] the solid line is the combined univer- 
sal law and the recommended super-catastrophic power 
law (equations [5] and [44| . The dotted line is our fit to 
disruption data on solid ice from Kato et al.| (1992), 
M h /M tot = 0.125(Q fl /QV 145 



The dashed line 



is our fit to disruption data on basalt from Fujiwara] 
et al. 1 1,1 °' 7 ' 7 ' 1 fl/f - /A/r — n A K -' 7tn - /r >'* 
that 



(1977], Mi r /Mt ot = QAbl{Q R /Q% D )-^\ Note 
ab data are available up to very high values of 
Qr/Q'rd ~ 100. The lab data spanning very weak to 
very strong geologic materials can be considered lower 
and upper bounds for the parameters in equation |44| 

The general agreement between the gravity and 
strength regimes suggests that gravitational reaccumu- 
lation of fragments has a negligible effect in the super- 
catastrophic regime. In other words, the mass of the 
largest fragment is primarily controlled by the shatter- 
ing process. 

Based on the similarity of the size distribution of frag- 
ments in laboratory experiments to the gravity regime 
data presented here (Figure [5]), we suggest that the dy- 
namical properties of the smaller fragments in super- 
catastrophic collisions are similar to the disruption 
regime. Therefore, the size and velocity distributions 
described in §3.2|and §3.3|can be applied. 



4.2. Hit-and-run regime 

Non-grazing impacts in the gravity regime transition 
from perfect merging to the disruption regime with in- 
creasing impact velocity. However, for impact angles 
greater than a critical v alue, an intermediate outcome 



may occur: hit-and-run ( Agnor fc Asphaug 



phaug et al.l|2006l llVlarcus et al.||2009l \'Mdb[ lAsphaug 



2004al As 



2010; Leinhardt et al. 20101. In a hit-and-run collision 
the projectile hits the target at an oblique angle but sep- 
arates again, leaving the target almost intact. Some ma- 
terial from the topmost layers of the two bodies may be 
transferred or dispersed. Depending on the exact im- 
pact conditions, the projectile may escape largely intact 
or may sust ain significant d amage and deformation (e.g., 
Figure 7 in |Asphaug||2010[ ). 
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Fig. 9. — Accretion efficiency (equation |45| versus velocity at infinity normalized by mutual escape velocity for different projectile-to- 
target mass ratios and impact parameters. Note that the impact velocity V% = y V^ f + V^ BC . Results from this work are connected by solid 

lines; previous results for supersonic impacts between protoplanets are connected by dashed lines l |Agnor fc Asp haug 2004 a|b| l and symbols 
are an aid to differentiate simulation groups. Magenta lines are for b = 0.5. 

The hit-and-run regime is clearly identified by consid- 
ering the accre tion efficiency of a collision, defined by 
Asphaugl 120091) as 



parameters. Data for colli sions between protoplanets at 
supersonic vel ocities from Agnor fc Asphaug| ( |2004a|b| ) 



Mr - M 



targ 



M n 



In a perfect hit-and-run event (Mi r = M targ ), £ = 0. 



M 



targ 



M p ), ? 



(45) 

For 
= 1. 



a perfect accretion event (Mi 

An erosive event in which Mi r < M tar g leads to £ < 
Note that the negative value of £ that corresponds to 
catastrophic disruption (Mi r = 0.5M tot ) depends on the 
specific mass ratio of the two bodies (£* = 0.5 — O.5/7). 

There is remarkably good agreement in the accretion 
efficiency and transitions from merging to hit-and-run 
and from hit-and-run to disruption between this work 
and previous simulations of higher velocity impacts be- 
tween large planetary bodi es ( Agnor_fc Asp haug 2 004a|b 



Marcus et al.|2009 2010b I. Figure [9] shows the accretion 



efficiency from our simulations in solid colored lines for 
four different projectile-to-target mass ratios and impact 



(and plotted in |Asphaug|2009| are shown in dashed lines 
for the common mass ratios (1:1 and 1:10). Hit-and-run 
collisions are indicated by a sudden drop from merging 
outcomes (£ = 1) to a nearly constant value of £ ~ for 
a range of impact velocities. Note that the drop in £ is 
sharpest for equal-mass bodies. For smaller mass ratios, 
the transition is not as sharp, and partial accretion of the 
projectile occurs at energies just above perfect merging 
(£ just above 0). 

Outcomes that are defined by the disruption regime 
have steep negative sloped accretion efficiencies. The 
disruption regime equations apply for partial accretion 
(0 < £ < 1) and for erosion of the target (£ < 0). Note 
that for high impact parameters (e.g., b = 0.7), there ex- 
ists an intermediate regime where the accretion efficiency 
has a very shallow negative slope and values of £ just be- 
low 0. These impact events, termed erosive hit-and-run, 
lead to some erosion of the target and more severe defor- 
mation of the projectile. The erosive hit-and-run regime 



1G 



Leinhardt and Stewart 



is eventually followed by a disruptive style erosive regime 
at sufficiently high impact velocities. The post-hit-and- 
run disruptive regime may be identified by finding the 
impact energy that leads to a linear relationship that 
satisfies the universal law. The required impact velocity 
increases substa ntially with in c reasing impact parame- 
ter; see £[5] and Marcus et al. (2010b) for an example 
disruption regime after a hit-and run regime (7 = 0.5 
and b = 0.5). 

In an ideal hit-and-run event, the target is almost un- 
affected by the collision, and the velocity of the largest 
remnant (the target) is about equal to the initial speed 
of the target with respect to the center of mass. More 
commonly, there is a small velocity change in both bod- 
ies which increase s the probability of merging in subse- 
quent encount ers ( Kokubo fc Genda|2010 ). Agnor & As- 
phaug ( 2004a]) referred to this collision outcome as inelas- 
tic bouncing. In our hit-and-run simulations with b = 0.9 
(green cluster of points in Figure [6] at M\ T = Mt arg ), the 
targets typically lose about 10% 01 their pre-impact ve- 
locity. For b = 0.7 (red points), there is more significant 
slowing of the target. Our data does not provide a robust 
description of the dependence of the post-impact velocity 
on the impact parameter and impact velocity. 

The projectile may be significantly deformed and dis- 
rupted during a hit-and-run event. The level of disrup- 
tion of the projectile may be approximated by consider- 
ing the reverse impact scenario: a fraction of the larger 
body impacts the smaller body. In this case, we esti- 
mate the interacting mass from the larger body with a 
simple geometric approximation. For the example geom- 
etry given in Figure [2j the cross-sectional area of the cir- 
cular projectile interacting with the target is calculated. 
The apothem is given by I — r, and the central angle 
is 4> — 2cos~ 1 ((Z — r)/r). Then, the projectile collision 
cross section is 

^interact = ?" 2 (7r - (0 - SU10)/2) (46) 

The interacting length through the target is approxi- 
mated by the chord at 1/2, 

L interact = 2 v AR 2 -(i?-V2) 2 - (47) 
And the interacting mass from the target is of order 



The reverse variables are 



.4 



interact -^interact 



(48) 



Note that the interacting mass depends on the impact 
angle (through I). 

To estimate the disruption of the projectile, we con- 
sider an idealized hit-and-run scenario between gravity- 
dominated bodies: the fraction of the target that does 
not intersect the projectile is sheared off with negligi- 
ble change in momentum and gravitationally escapes the 
interacting mass. Hence, we ignore the escaping target 
mass and consider only the impact between Mi n t orac t and 
the projectile mass, M p . 
The reverse impact is thus defined by = Mi n t orac t 

t denotes the reverse im- 
For each impact angle, calculate R} C1 
from the principal dis- 



and A/ t T arg = M p , and the 



pact variables. 



for M, 



Aft 



mL s , Q 



t* 

_RD,7=1 



ruption curve (equation 28), and Vl^ from equation 



30 



MlMlj{Ml 



7 t =Mt/M t t 



Ml 



(49) 
(50) 



'p/ ^targ- 

The mass ratio correction from the principal disruption 
curve is 

l/(3/Z) 



^RD 



1 ( 7 t + If 



7 



t* 



it* 



l(7t 
4 



7' 



2/(3/2)-l 



(51) 



(52) 



Once the reverse impact disruption criteria is calculated, 
we use the universal law for the mass of the largest rem- 
nant to determine the collision regime for the projectile. 
If the projectile disrupts, then the size distribution of the 
projectile fragments may be estimated to first order from 
the disruption regime scaling laws. 

5. TRANSITIONS BETWEEN COLLISION REGIMES 

5.1. Empirical transitions between accretion, erosion, 
and hit-and-run 

We have classified the collision outcome regime for all 
of our new simulations. The outcome is sensitive to the 
mass ratio of the two bodies, the impact parameter, and 
the impact velocity. Four regimes are mapped in Fig- 
ureHni 

1. Accretion of some or all of the projectile onto the 
target (M lr > M targ and £ > 1, light blue squares), 

2. Partial erosion of the target (M\ r < Af targ and £ < 
1, dark blue squares), 



Mt ar g and £ = 0, green 



3. Pure hit-and-run (Mi r 
triangles), and 

4. Erosive hit-and-run (M\ r slightly less than M targ 
and £ slightly less than 0, red triangles). 

Note that the 1:40 mass ratio simulations reach im- 
pact velocities that exceed the physics included in PKD- 
GRAV; impact velocities greater than about 1 km s _1 
should use a shock physics code. Hence, the transition 
to the erosive regime at high impact parameters could 
not be derived directly. 

For impacts at small impact parameters (more head- 
on), the collision outcomes transition from accretion to 
erosion with increasing impact velocity. For more oblique 
impacts, the collision outcomes transition from merging 
to hit-and-run to eros ion with increas ing impact veloc- 
ity. As suggested by Asphaug (2010), 6 C rit (red verti- 
cal line in Figure 10 ) is indeed a good indicator of the 
minimum impact parameter necessary to enter the hit- 
and-run regime. However, for 7 = 1, we find a small 
region of erosive hit-and-run events when b < b cl -a . The 
use of 6 cr it to define grazing and non-grazing impacts 
makes the very simplifying assumption that the velocity 
vector of the center of mass of the projectile remains con- 
stant during the event. In reality, the projectile center 
of mass will be deflected to some extent during the en- 
counter, and the true interactive mass will be larger than 
assumed here. The deflection is greatest for more equal- 
mass bodies, and a narrow region of erosive hit-and-run 
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Fig. 10. — Map of the major collision regimes as a function of mass ratio, impact parameter, and impact velocity normalized by the escape 
velocity from the combined mass with radius Rqi- Cyan squares — a full or partial accretion event, M\ x > Mtarg; blue squares — target 
is eroded, M\ T < Mtarg; green triangles — ideal hit-and-run event, M\ T = Mtarg; red triangles — erosive hit-and-run event, M\ r slightly 
less than Mtarg- Red vertical line co rresp onds to 6 cr it for the given mass ratio (equation Black curve is onset of erosion predicted from 
the catastrophic disru ptio n model (j ]3.1| with c* = 4.3 and p, = 0.35; dashed black curve is predicted transition from perfect merging to 
hit-and-run (equation |53|> - 



events is observed for b = 0.35 and 7 = 1. Note that the 
transition between erosion and hit-and-run occurs near 
b CT it for all size bo dies studied to date, fr om 1 km rubble 
pile planetesim als ( Leinhardt et al.|2000 ) to super-earths 
( |Marcus et aL|2009| r~ 

For grazing collisions, the hit-and-run regime is 
bounded by perfect merging at low impact velocities the 
onset of disruption at high velocities. The projectile 
merges with the target when the impact velocity is less 
than the mutual escape velocity (in other words, the ve- 
locity at infinity Vinf is zero). Since only a fraction of the 
projectile may interact in oblique impacts, the appropri- 
ate mutual escape velocity for perfect merging is slightly 
less than the mutual escape velocity from the total mass. 
Then the appropriate measure for merging is 



K'sc = y/QGM'/R'), 



(53) 



where M' = 
((3M')/(4tt /9 )) 1 /3 i 



targ 



T" ^-interact 



and B! 



target have the same bulk density p. The boundary 
between me rging an d hi t-and-run is well matched by 
equation 



53 in Figure 10 (dashed black line). 



assuming that the projectile and 



Of course, the concept of an interacting mass is a sim- 
plistic limit because it assumes that the part of the pro- 
jectile that impacts the target can separate from the rest 
of the projectile without loss of momentum. In Figure [T0| 
the only set of simulations that did not show a sharp 
transition from merging to hit-and-run is 7 = 0.1 and 
b = 0.7. In this case, the impact parameter is very 
close to 6 cr it = 0.66, and the outcomes include partial 
accretion of the projectile, erosive hit-and-run, and fully 
erosive collisions with increasing impact velocity. 

Loss of momentum by the projectile in grazing colli- 
sions does lead to merging when Vi n f > 0; in Figure [9j 
note the nearly complete merging in the 7 = 1 sim- 
ulations fo r small values of Vj n f with b = 0.5 and 0.9 
(data from Agnor fc Asphaug| |2004a |) . For Vi n f slightly 
abov e zero, merging occur s m graze-and-merge events 
(e.g., Leinhardt et al.pblO ). In such cases, the two bod- 
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ies hit, separate as nearly intact bodies with decreased 
velocity, and then merge upon a second collision. The 
impact velocity range for graze-and-merge outcomes is 
quite narrow; previous studies have demonstrated the 
small velocity increase needed to transition from perfect 
merg i ng to graze-and-merg e to hit-and-run (e.g., Canup 



2004 



Leinhardt et al.|2"0 10). Concurrent with this work, 
the graze-and-merge regime has been explored i n more 
detail using hydrodynamic SPH simulations by |Genda 



et al.| ( |2011b|) 

Grazing collisions transition out of hit-and-run to ero- 
sion of the target when the impact velocities reach the 
disruption regime. The transition to the disruption 
regime is a strong function of the impact parameter be- 
cause of the rapidly shrinking projectile interaction mass 
and the dependence of the disruption criteria on the mass 
ratio and impact velocity. 

Our general model for the catastrophic disruption cri- 
teria combined with the universal law for the mass of 
the largest remnant is used to derive the impact velocity 
needed to begin ero ding the target mass (M\ T — M tar g, 
black line in Figure 10). Our new simulation data are 
best fit with a value oFTr = 4.3 and p, = 0.35. Our model 
for the disruption regime provides an excellent estimate 
for the transition to erosion of the target for the wide 
range of impact parameters considered here. In particu- 
lar, the analytic model captures the sharp increase in the 
upper bound to the hit-and-run regime between b = 0.7 
(45°) and 0.9 (64°). 

5.2. Predicted transitions between accretion, erosion, 
and hit-and-run 

Using our analytic model, we derive example collision 
outcome maps for collisions between protoplanets. We 
fit values of p, = 0.36 and c* = 1.9 to the data from 
collisions between planet-sized bodies using SPH codes 
(Figure |4^3). Collision maps, which are color-coded for 
outcome regime, are shown in Figure 11 for four mass 
ratios. 

The details of the forward calculation of the collision 
regimes are given in the AppendijQ In Figure 11 the 
impact parameter axis is scaled by the probability of 
an impact at that angle. The probability of an im- 
pact within a n interval (9 , 9 + d 9) is proportional to 
sin(0) cos(f?)df? Shoemaker (1962). The corresponding 
impact angle is shown on the top axis with 5° tick inter- 
vals. The model assumes an abrupt transition between 
grazing and non-grazing impacts, which is certainly arti- 
ficial. Near the critical impact parameter, collision out- 
comes will have elements from both the disruption and 
hit-and-run regimes. 

Contours of impact velocities that correspond to a con- 
stant mass of the largest remnant are calculated using the 
general model for catastrophic disruption and the univer- 
sal law or powe r la w for the m ass of the largest remnant 
(equations [5] or 44). In Figure 11 the thick black curve 
corresponds to the critical velocity for catastrophic dis- 
ruption, where the largest remnant contains half the to- 
tal mass. Note that this curve corresponds to the target 
erosion boundary for 1:1 scenarios (the transition from 
partial accretion (light blue) or hit-and-run (green) to 

4 A code to generate collision outcome maps and to calculate 
specific impact scenarios is available from the authors. 



erosion (white) regions). The grey dashed curves corre- 
spond to the impact velocity needed to disperse 10% and 
90% of the target mass. 

Between perfect merging and erosion of the target, 
there is a region of partial accretion of the projectile. For 
non-grazing impacts, the dotted curve corresponds to ac- 
cretion of 50% of the projectile mass. Grazing impacts 
transition rapidly between perfect merging and hit-and- 
run with increasing impact velocity. 

Most hit-and-run collisions with M p < O.lMtarg result 
in significant disruption of the projectile. In the collision 
outcome maps, the onset of projectile erosion in a hit- 
and-run event is given by 



(54) 



targ 



Note that for impact parameters near & cr j t , M_[ ~ M i 
for projectile-to-target mass ratios less than about 0.1. 
Thus the velocity contours of constant remnant mass in- 



tersect for catastrophic disruption (MJ = 0.5M t T ot ) and 
onset of projectile erosion (M^ = M targ ). Futhermore, 
there is a minima in the projectile erosion curve at an 
optimal fraction of total interacting mass from the tar- 
get (in other words, the reverse projectile-to-target mass 
ratio is varying with impact parameter). The two veloc- 
ity contours diverge at higher impact parameter as 

becomes much less than M t ' arg . 

Coll ision maps for planetesimals are presented in 
j 6.2.1 and the implications of the diversity of colli sion 
outcomes for planet formation are discussed in §6.3[ 

6. DISCUSSION 

6.1. Scaling of collision outcomes in the gravity regime 

For all gravity-regime bodies studied to date, collision 
outcomes are strikingly similar for a tremendous range 
of target composition and size. Furthermore, the transi- 
tions between the major collision regimes (merging, hit- 
and-run, disruption, and super-catastrophic disruption) 
occurs under similarly scaled conditions. The types of 
bodies studied, ranging from km to several 1000's km in 
size, i ncluded rubble-pile and porous planetesimals (this 
work. 



2010 



Stewart fe Leinhardt||2009l |Benz||2000| | Jutzi et al. 
Korycansky & Asphaug 2009), pure rock or pure 



ice pla netesimals with stre ngth ptenz fc Asphaug 1999 
Leinhardt fc Stewart||2009l [jutzi et al.||2010|), sTrengflT- 



le 



ess differentiated rock and iron plane ts" jfjenz et al. 
19881 [20071 1 Agnor fe Asphaug|2004a| |Marcus et, alJiMW 
Gend a et al.|2011b[), strengthless d ifferentiated water ani 



rock planets (|Marcus et al.|2010b|) , and strengthless pure 



rock planets ( |Marcus et al. 2009p . The studies focused 
on a variety of stages during planet formation, from ac- 
cretion of planetesimals to destruction of planets; thus, 
the impact velocities ranged from ~ 1 m s _1 to over 100 
km s . The computational methods included three dif- 
ferent shock physics codes and two TV-body codes. Our 
analysis of the results from these studies suggest that 
the same scaling laws may be applied over an incredibly 
broad range of im pact scenarios d uring planet formation. 

As stressed by Asphaug (20101, similarity of outcome 
is not the same as true scale invariance. He notes that 
scale invariance applies only for idealized incompressible, 
self-gravitating inviscid fluid planets. In reality, many 
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Fig. 11. — Predicted collision outcome maps using the analytic model for strengthless planets (c* = 1.9 and fi = 0.36) for selected 
projectile-to-target mass ratios. Impact velocity is normalized by the mutual surface escape velocity assuming a bulk density of 3000 kg m — 3 ; 
impact parameter is spaced according to equal probability. Colored regions denote perfect merging (dark blue), partial accretion (light 
blue), net erosion to the target (white) and hit-and-run (green). Vertical red line denotes the onset of hit-and-run events at 6 cr it. Thick 
black curve - critical disruption velocity for half the total mass remaining; grey dashed curves - 10% and 90% of target mass in largest 
remnant; dotted curve — 50% of projectile accreted; dot-dashed blue curve - catastroph ic disruption of the pro jectile; dashed blue curve — 
erosion of the projec tile. Example pr oposed giant impact events: • — Haumea system l |Leinhardt et al.|2010[ |; A - Pluto-Charon | |Camip| 
2005); ♦ - Mercury | |Benz et al.|2007| l; ■ - Earth-Moon JCanup|2004) . 



aspects of collision outcomes will not scale simply with 
size: e.g., the mass of collision-produced melt depends 
on the specifics of impact velocity, target composition, 
and the internal temperature and pressure history. Here, 
we investigated the similarity of the dynamics of colli- 
sion outcomes for a variety of non-ideal gravity-regime 
bodies, from icy planetesimals to differentiated super- 
earths. Specifically, we developed scaling laws to define 
the mass and velocity distribution of bodies after any 
gravity-regime collision. 

Why do the dynamics of collision outcomes appear to 
scale similarly with size in the gravity regime? At im- 
pact velocities just above the escape velocity, momentum 
dominates the outcome at all scales. Hence, the transi- 
tion from merging to hit-and-run depends primarily on 
the geometric cross section of the collision for all size bod- 
ies. As impact velocities increase, the energy required for 
disruption is dominated by the gravitational dispersal of 
fragments rather than the en ergy required to shatte r an 
intact body into small pieces ( Melosh fc Ryan|1997[ ). As 
a result, erosive outcomes require that the velocity of 
the fragments exceed a critical value that relies primar- 
ily upon the gravitational potential of the total colliding 
mass. 



For small bodies, the critical fragment velocity may be 
reached with impact velocities that impart negligible ir- 
reversible work on the materials (Figure |4]) . For larger 
bodies, the critical velocity requires sufficiently high im- 
pact velocities that strong shock waves are formed. The 
shock wave permanently deforms the materials and, in 
the process, reduces the total energy available for the 
final velocity distribution of fragments. The energy of 
deformation is often referred to as "waste heat"; for a 
fixed impact energy, a larger fraction of waste heat is 
generated with increasing impact velocity (primarily due 
to the onset of shock-induced melting and vaporization 
at high shock pressures). As a result, the catastrophic 
disruption criteria increases with increasing impact ve- 
locity (equation [3]) . 

Based on currently available data we argue that in the 
disruption regime the dynamics of the outcome is sim- 
ilar over the entire gravity regime when scaled by the 
catastrophic disruption criteria. The post-collision size 
distribution is similar, as it is controlled by the largest 
remnant and the gravitationally accreted clumps from 
the shattered parent bodies. The general catastrophic 
disruption law accounts for both the increasing gravita- 
tional potential with total mass of the colliding bodies 
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and the increase in waste heat at higher impact veloci- 
ties ( |Housen fc Holsapple||1990[ ). 

The development ot equation [3] relied upon the con- 
cept of a coupling parameter, A oc RpV/ 1 , a point source 
approximation of the coupling of t he projectile's energy 
and m omentum into the target (c.f.lHolsapple & Schmidt 
19871. The velocity exponent p is bounded by pure mo- 
mentum coupling (p = 1/3) and pure energy coupling 
(p = 2/3). In the gravity regime, the coupling parame- 
ter distills the physical response of the geologic material 
into the variable p. Some constraints on p are avail- 
able from laboratory cratering exp eriments: e.g., p = 0.4 
for sa nd and p = 0.55 for water (Holsapple & Schmidt 



19871. Here, we fit the coupling parameter to the numer- 
ical simulation results for disruption of a wide variety of 
materials. The derived best fit range of 0.33 < p < 0.37 
is close to pure momentum scaling. 

Why does the concept of a point source approximation 
apply to collisions between comparably sized planetary 
bodies? The point source approximation was developed 
for impact crat ering by a finite size pr o jectile onto a half 



space target. Holsapple & Schmidt (1987) show that 



the concept of a point source is equivalent to a variety 
of models that describe a similar material velocity field 
far from the impact point. In the case of catastrophic 
disruption, the late-time far-field criteria is a fragment 
size- velocity distribution where half the mass is escaping 
the gravitational potential of the largest remnant. The 
principal dynamical factors governing the collision are 
incorporated into the Qr D formulation: relative veloc- 
ity, mass ratio, impact parameter, and bulk density. The 
similar outcomes of collisions with similar Qr/Q'rd in- 
dicate that the remaining details of how the energy and 
momentum is distributed into the target and projectile 
during the initial stage of the collision are negligible in 
determining the late time dynamics following a catas- 
trophic disruption event. 

In summary, the primary factors that bound the differ- 
ent collision outcomes regimes scale similarly with size in 
the gravity regime: momentum, geometric cross section, 
and normalized critical impact energy (Q'rd)- Other fac- 
tors that lead to second order perturbations to the dy- 
namics of the collision outcomes are discussed in the next 
section. 

6.2. Scaling laws limits of applicability 

Planet formation involves a vast range of bodies with 
distinct physical characteristics, including dust aggre- 
gates, rubble-pile planetesimals, differentiated molten 
and solid protoplanets, solid planets with extended at- 
mospheres, and gas-dominated planets. The constituent 
materials (iron-alloys, silicates, ices, and gases) span or- 
ders of magnitude in density and material strength. The 
complex and time-varying physical properties of plan- 
etary bodies significantly limits the application of any 
single equation to all bodies over the course of planet 
formation. And yet, judicious simplification is necessary 
for planet formation simulations to be both physically 
robust and computationally tractable. 

We have focused on developing scaling laws that de- 
scribe the dynamical outcome of collisions between any 
two gravity-dominated bodies. The dynamical outcome 
from collisions seems to be rather insensitive to the in- 
ternal composition, when the results are scaled by the 



appropriate value for Q* RD - However, the types of bodies 
studied to da te do not contain any g as mass fraction (see 
for example, Kobayashi et al. 2011), and so the scaling 
laws may need modification for a planet with a signifi- 
cantly different internal structure than the differentiated 
and homogeneous planets included in this study. One 
area that warrants further investigation is the sensitiv- 
ity of the velocity of the largest remnant to the internal 
structure / composition (Figure [6| . 

The role of tidal effects during collisions or in close en- 
counters may be importa nt factors during the fragmen- 
tation of planetary bodies ( Asphaug et al.|2006j|Asphaug| 
2010 1. In this work, all bodies are assumed to be approx- 



imately spherical at the time of impact. Tidal affects will 
change the interacting mass and contribute to the frag- 
mentation process in ways that lead to different size and 
velocity distributions than found here. 

Similarly, the role of pre-impa ct spin during collisions 



has received modest attention (Leinhardt et al. 2000 



Canup|2008 1 . While the net spin of a growing body may 



essentially average near zero during the rapid growth 
phase where collisions are frequent, the effect of pre- 
impact spin and the collision angular momentum may be 
very impor tant in the final gia nt impact phase of planet 
formation (Agnor et al. 19991. While, we did not con- 
sider any pre-impact spin in this study, a few simulations 
with high collision angular momentum are notable. In 
1 : 1 collisions w ith < b < 6 cr jt (Tab le Jl] in this work 
and Table 1 in Leinhardt et al. |2000[) , the catastrophic 
disruption criteria is less than the value at b — (e.g., 
closed and open stars in Figure |4| . We interpret the 
lower disruption criteria with pre-impact spin to arise 
from the significant collision angular momentum. As a 
result, the gravitational potential is reduced and disper- 
sal requires slightly less energy. We suggest that future 
work investigate the possibility of using the reduction in 
the gravitational potential due to pre-impact and colli- 
sion angular momentum to account for the first order af- 
fects of spin. Specifically, the spin-modified catastrophic 
disruption criteria may remain at a constant offset (c*) 
from the spin-modified gravitational potential. 

6.2.1. Strength and porosity in the gravity regime 

The study of catastrophic disruption of strong rock tar- 
gets has been motivated by collisional evolution studies of 
the asteroid and Kuiper belts. The strength models were 
tested by fitting laboratory quasi-static strength mea- 
surements and fragment size distributions from head-on 
disruption experiments. Particular attention was p aid to 
the development of the model for tensile fracture (Benz 
& Asphaug 19941, as the tensile strength dominates the 
catastrophic disruption criteria for head-on impacts in 
the strength regime. 

Results from several numerical simulations of catas- 
trophic disruption of strong rock targets in the grav- 
ity regime are shown in Figure [12] The head-on basalt 
disruption data at impact velocities of 3 a nd 5 km s -1 
(*) are show n from the canonical study by Benz & As 



phaug ( 1999 ) using the SPH code with the detail e d ten- 
sile strength model. Using the same code, Benz (2000) 
studied the disruption of strong nonporous basalt (hour- 
glass) and a macroporous target, composed of overlap- 
ping clusters of SPH particles representing strong inter- 
connected boulders (m), at very low impact velocities 
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Fig. 12. — Catastrophic disruption simulation results for strong 
rock targe ts (p orous and nonporous). Same notation as in FigureRl 
and Table |AT2] A. Simulation data corrected to an equivalent heao> 
on impact! B. Simulation data converted to an equivalent equal- 
mass disruption criteria. The results for critical velocities from 
m s _1 to 5 km s — 1 demonstrate that energy scaling is incorrect. 
Best fit Q* RD curves with p. = 0.35 and c* = 6.4 for V* = .005, 
.02, .1, .3, 1.5, and 5 km s" 1 . 



(5-40 m s _1 ) and b = 0.7. In Figure 12 the 10-km target 
data, which fall below the specific gravitational binding 
energy, are derived from the equal-mass collisions pre- 
sent ed in Benz's Figure 5 and will be discussed below. 
The Benz (2000) 1-km data are less certain using our 
catastrophic disruption variables because both impact 
velocity and mass ratio were varied and the specific val- 
ues were not reported. Nevertheless, the significant offset 
in the disruption criteria is irrefutable evidence that pure 
energy scaling does not apply. In fact, the total disper- 
sion in the specific impact energy is slightly larger than 
can be accommodated by the momentum scaling limit 
of p = 1/3, which is likely a result of differences in the 
details of the strength models. 
More recent simulations (<£). 



with critical velocities of 3 and 5 km s 



by |Jutzi et al.| ( |2010 | 
1 fall in-between 



the data from |Benz fc Asphaug (1999). Their work uses 
the same SPH code with an updated strength model that 
includes the extra dissipation of microporosity. Simula- 
tions using the CTH shock physics code with different 
shear and tensile strength models yield similar results as 
fou nd for nonporous basalt tar gets using the SPH code 
(▼, |Leinhardt_fc Stewart|2009| ). 

In Figure |12| the strong target data, with mass ratios 
from 1:1 to almost 1:45,000 and impact velocities from 
0.001 to 5 km s" 1 , are best fit by p = 0.35 and c* = 6.4. 
The equivalent equal-mass disruption data have c* val- 
ues from 1 to 20. For comparison, the best fit to only the 
PKDGRAV rubble pile data is c* = 5.5 and p = 0.365. 



We note that the data from Jutzi et al. (2010) and the 
1-km targets from Benz (2000) nic ely cluster ar ound the 
best fit principal disruption curve. Jutzi et al. (20101 fit 
their 3 and 5-km s" 1 data with p = 0.43; however, such 
a high value for p cannot simultaneously fit the data at 
lower velocities. The two-dimensional simulations from 



Leinhardt & Stewart ( 2009 ) fall systematically below the 
best fit curve. The 3 an d 5-km s _1 head-on data from 
Benz fc Asphaug] ( |1999[ ) have a dispersion greater than 
can be explained with our model; the low and high- 
velocity data fall below a nd above the best fit cu rve, 
respectively. The da ta from Benz fc Asph aug (1999) and 
the 10-km data from [Be nz (20 00) were excluded from the 
global fits presented in J 3.1.3' 



The 10-km equal-mass data from Benz (2000[) (ix and 



hourglass) require closer examination. At an impact an- 
gle of 45 degrees, all 1:1 data on weaker bodies pass 
through the hit-and-run regime. However, both the non- 
porous and porous data show disruption results similar 
to the non-grazing regime. We interpret the non-grazing 
outcome to be due to the high shear strength of the tar- 
get preventing a hit-and-run outcome. We hypothesize 
that the disruptive outcome and disruption energy be- 
low the gravitational binding energy are related to the 
strength and angular momentum of the event. A colli- 
sion between two equal-size strong bodies has a larger 
interacting mass than assumed in our model, so the ad- 
justment from the oblique to equivalent head-on colli- 
sion disruption energy is overestimated. In addition, the 
collision generates significant spin angular momentum. 
The angular momentum reduces the effective gravitional 
binding energy and, similarly, the required disruption en- 
ergy. These data illustrate the need to better understand 
the physical properties of strong targets in oblique im- 
pacts and the role of angular momentum. 

In the strong rock target simulations, the typical lim- 
iting shear strength is 3.5 GPa, comparable to the quasi- 
static shear strength in laboratory rock under high con- 
fining pressure. In the SPH simulations, the shear 
strength was fixed throughout the simulation. In the 
CTH simulations, the shear strength was dependent on 
the confining pres sure and the accumula t ion of damage 
(e.g., fractures). Leinhardt fc Stewart| (2009) demon- 
strated that shear strength is important for the size bod- 
ies considered here, which are usually considered to be 
purely in the gravity regime. Higher shear strength leads 
to greater dissipation of the shock energy into material 
deformation; hence, higher specific energies are required 
to disrupt stronger targets. None of the published work 
has investigated the role of strain rate on zones of shear 
localization in catastrophic disruption simulations, which 
leads to significant reduction of shear strength durin g 



impact cratering events (e.g., Senft fc Stewart||2009[ ). 
More work is needed to develop more sophisticated shear 
strength models for strong rock targets and to validate 
model calculations for oblique impacts. 

There has been some recent work on the catastrophic 
disruption of porous planetesimals. Porosity has been 
modeled in three different ways: as hard sphere rub- 
ble piles with various bulk densities in studies using 
PKDGRAV (see references in Figure [2]), macroporous 
overlappi ng clusters o f SPH particles representing intact 



boulders (Benz 2000), and microporous bodies using a 
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constitutive model for porosity in an SPH code (Jutzi 
et al. 20101. The SPH simulations found significant ef- 
fects of porosity in the strength regime; however, poros- 
ity was a second order effect in the gravity regime, and 
the catastrophic disruption criteria agreed with the non- 
porous simulations when the data were normali zed by 
the difference in bulk density (Jutzi et al. 2010). The 
low-velocity macroporous SPH simulation results in the 
gravity regime o verlap with the PK DGRAV rubble pile 
results. Finally, |Jutzi et al.| (2010) found similar frag- 
ment size and velocity distributions between their porous 
and nonporous gravity regime results. 

We note that the transition between the gravity and 
the strength regime should be handled carefully and ap- 
propriate coefficients should be chosen for different mate- 
rial composition and strength. There appears to be sig- 
nificantly more variation in the disruption criteria in the 
strength regime compared to the gravity regime; how- 
ever, future work should consider whether or not a scal- 
ing analysis similar to the one presented here may cap- 
ture most of the variance. 

6.2.2. Other collision outcomes 

In cases where the impact velocity is above the escape 
velocity but the mass of the projectile is too small to lead 
to disruption, some material will esca pe the target in the 
form of cr ater ejecta. In recent work, |Housen fe HolsajT] 
pie (2011) has conducted a detailed study of the scaling 
of ejecta from im pact craters. Based on man y labora- 
tory experiments, |Housen fc Holsapple (2011) find that 
approximately 0.01iW p of material achieves escape veloc- 
ity in cratering events at Vi ~ V^ sc (see their Figure 16). 
Empirical fits to the material eroded during cratering 
events onto sel f-gravitating b odies has also been studied 
numerically by Svetov (2011). 

The bulk composition of a, body may change during 
planet formation by either preferentially accreting mate- 
rial of a certain composition (e.g., ice fragments chipped 
off smaller bodies) or by stripping of mantle material. 
The loss of a mantle during catastrophic disruption has 
been studie d for rock/iron and w a ter /roc k differentiated 
planets by |Marcus et al] ( |2009| |2010b[ ). They devel- 
oped two models to calculate the resulting change in the 
mantle mass fraction that bound the simulation results. 
Their method for calculating the change in composition 
is described in the Appendix and may be incorporated 
into planet formation studies that track the composition 
of gr owing and eroding planets (Stewart & Leinhardt 
20111. 

Previous work has addressed collision outcomes in the 
strength regi me to various levels of gene ralit y. We refer 
the reader to|Beauge fc Aarseth (1990) and Kenyon & 
Bromley ( 20081) and references therein. 



b odies (|AlphauggJal.j2006| |Asphaug||2010[ ). 
' Agnor~et al.| ( 1999 ) found that the most common col- 



6.3. Implications for planet formation 
6.3.1. Giant impact events 

Even with limited understanding of the full dynam- 
ics of collision outcomes, the significant role of giant im- 
pact events in determining the final physical p roperties of 
rocky /icy planets has been recognized (e.g., Agnor et al. 
1999| Ida fc Lin||2010 ). Any event between similar sized 
bodies (7 > 0.1) may be considered a giant impact event, 
although the outcomes are more dramatic for larger mass 



lision events at the end stage of terrestrial planet for- 
mation (under our solar system conditions) have 7 ~ 
0.01 — 0.2 and Vi from about 1 to 4V r 0SC . Over this range 
of mass ratios and impact velocities, collision outcomes 
span all the re gim es: accretion, erosion, and hit-and- 
run. In Figure nip, note that 6 cr i t = 0.66 for 7 = 0.1; 
hence, about half of all impacts fall in the regime that 
transitions from accretion to erosion and half transition 
through a hit-and-run regime. Hence, the implementa- 
tion of self-consistent scaling laws to describe collision 
outcomes is crucial to the accurate treatment of the gi- 
ant impact phase of planet formation. Although Agnor 



et al. ( 1999 ) typically found that impact velocities fell in 
the range of 1 to 2V CSC , temporary dynamical excitation 
by migrating giant planets may further increase the im- 
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robust character] 


zation of all 



collision outcomes is necessary for any planet formation 
calculation. 

With the strong dependence of collision outcome on the 
mass ratio, the final stage of planet formation is likely 
to produce more diverse outcomes than previously an- 



ticipated. As argued by Asphaug (20101, the increased 
frequency of hit-and-run events during the giant impact 
stage may routinely lead to compositional modi fica tion 
of the second-largest body. As shown in Figure [TT] the 
escaping projectile is nearly always eroded in hit-and- 
run events. Consequently, the atmosphere, hydrosphere, 
and even the mantle of these bodies may be stripped 
away. Such interesting details may no w be explo r ed di- 
rectly in planet formation simulations. lAsphaug (20101 



suggested that the growth of large rocky planets occurs 
often by a series of hit-and-run events followed by an 
eventual merger. Under these circumstances, each ac- 
creting protoplanet could have been partly devolatilized 
before merging. In this manner, the final composition of 
planets may be altered significantly compared to the ini- 
tial protoplanets during accretion into the final planets. 

Note that our analytic calculation of collision outcomes 
agrees very well with the proposed giant i mpact scenarios 
for the formation of the Haumea system ( |Leinh ardt et al. 
2010 1, the formation of Pluto a nd Charon (|C anup|2005 ) , 



the formation of Earth's moon (|Canup | 2004 |), and the in- 
creased density of Mercury ( |Benz et al.|2007p (Figure 11 ). 
The formation models plotted tor Haumea and Pluto- 
Charon are the result of graze-and-merge events, where 
two equal-mass bodies collide and separate nearly intact. 
The loss of velocity by the first collision leads to a merg- 
ing upon a second collision, creating a final body with 
enough angular momentum to spin off a disk of material. 
In contrast, the canonical formation of the Moon involves 
a collision where the projectile is disrupted upon the first 
impact. The analytic calculation for disruption of the 
projectile agrees very well with these moon-forming stud- 
ies. Because the giant impact phase of planet formation 
is dominated by collisions slightly above the mutual es- 
cape velocity, the probability scaled axis in Figure [TT] 
emphasizes the high likelihood that the giant impacts 
will be on the bo undary of the merging and h it-and-run 
regimes (see also Stewart & Leinhardt (2011)). 

Given the range ot impact velocities found by [Agnor 
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Fig. 13. — Predicted collision outcome maps for small projectiles 
(M p : Mtarg = 1 l 200) using the analytic model for (A) a nominal 
planetesimal with c* = 5 and (B) best fit to asteroid family for- 
mation simulations with c* = 20. Axe s, colors and line notations 
are the same as defined in Figure |11| Proposed astero id family- 



forming ev ents with M\ T 
|et al.|2006| |. 



0.1 - 03iAitarg: ▼ - Karin jNesvorny] 



et al. ( 1999 ) in the giant impact stage (up to about 6V^ SC ), 
stripping the mantle from Mercury by a catastrophic 
disruption event is reasonably probable. Recently, col- 
lision outcomes alone have been used to define a limit to 
the possible density of super-earth mass exoplanets (1- 
lOM®). Based on the criteria to s trip off mantle mate- 
rial during catastrophic disruption, Marcus et al. (2010a I 
used the extremely high impact velocities required to cfis- 
rupt 1 to 10 earth mass planets to place an empirical 
limit to the iron fraction of a planet that has suffered 
a single catastrophic impact event. The prediction is 
in very good agreement with observat ions of rocky exo- 
planets (e.g., Kepler 10b and 55 Cnc e, Winn et al |2011 



Batalha et al. 2011). Consideration should be given to 
the potential for stripping mantles off the planets by ero- 
sive hit-and-run events: e.g., the smaller projectile has 
its mantle stripped but it is never incorporated into a 
larger body. 

6.3.2. Collisional evolution of small body populations 

The asteroid and Kuiper belts contain a valuable 
record of the dynamics of the solar system. The colli- 
sional evolution and dynamical clearing of these reser- 
voirs of small bodies has been modeled extensively (e.g 



Nesvorny et al.|200 2; Bottk e et al.|2002||Morbidelli et al. 
2008||Kenyon et al.|2008[ ). Observations of asteroid fam- 



ilies, in particular, have been important in the study of 
planetary dynamics and impact processes. Asteroid fam- 
ilies and their associated dust bands are believed to have 



formed in recent catastrophic impact events ( [Nesvorny 
et al. 2003 ). Simulations of asteroid disruption have been 



compared to observations of the size and velocity distri- 
bution of asteroid fa milies to derive possib le impact sce- 
narios. For example, Nesvorny et al. ( 2006 ) simulated the 
formation and dynamical evolution of the Karin family. 
Using the same SPH co de and strength model used by 
Benz & Asphaug ( 1999[ ), their best fit scenarios for the 

impact onto a 15- km 
0.1-0.2M targ 



Karin family involve a 5 to 7 km s _1 
target with a mass ratio of 1:200 and M\ T 
(Figure 13). In order to match this combination of im- 



pact energy and largest remnant mass with our analytic 
model, a value of c* of approximately 20 is required, 
which is significantly higher than the best fit value of 
6.4 derived here for strong targets. Figure [13] presents 
collision outcome maps for the best fit for airsmall bod- 
ies (c* = 5) and the very strong bodies inferred from the 
asteroid-family formation simulations (c* ~ 20). Catas- 
trophic disruption begins at impact velocities of 107V r csc 
for the generic small body, whereas values about 2 times 
higher are requ i red fo r the strong targets simulated by 



Nesvorny et al.| ( |2006[ ). 

In addition to studying individual collisions, the colli- 
sional evolution of the asteroid belt size distribution has 
been studied in detail (e.g. 



Davis et al. 1979 Bottke 



et al.||2005| |Morbidelli et al.||2009| . Such studies seek to 
understand the relative weight or collisional versus dy- 
namical clearing of the belt and the initial size distribu- 
tion of planetesimals. From our discussion of the role 
of strength during the evolution of planetesimals from 
weak aggregates to protoplanets, it is clear that a single 
catastrophic disruption criteria cannot apply at all times. 

Also, the influence of mass ratio on the disruption crite- 
ria will be important if the first planetesimals were born 
big. Morbidelli et al. ( 2009 ) argue that the observed size- 
frequency distribution ot asteroids is best fit by an initial 
population of planetesimals that were 100's km in size. 
(Note this suggestion has not been fully accepted as a 
requirement for the observed the size distribution of as- 
teroids (see WeidenschillingpOll ).) Two processes have 
been proposed to form km to 100 k m-scale initial plan- 
etesimals: turbulent concen tration (Cuzzi et aD |2008| 
and streaming instabilities ( |Johansen et al. 200y[ )~ A 
mechanism to form km-scale planetesimals is attractive 
because it could help bypass the so-called meter-size bar- 
rier, the size particle that radially drifts in toward the sun 
fas ter than it can gr o w. In the collision evolution model 



by Morbidelli et al. 



threshold is the angle- averaged 5 km s 



(|2009|), the catastrophic disruption 
1 constant veloc - 
Jl999l. 



ity curve for strong basalt from Benz fc Asphaug| 
Note, however, that the proposed mechanisms 10 torm 
km-scale planetesimals would generate weak aggregates 
of small (e.g., cm to m-size) particles. These aggregates 
will be comparable to simulations of weaker materials. 



Morbidelli et al. (2009) considered a collisional evolution 
simulation using a catastrophic disruption criteria that 
was a factor of 8 lower than the basalt criteria. They note 
that the size-frequency distribution was not significantly 
affected by the larger amount of collisional grinding in 
the one test case; however, they could not form the lu- 
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nar to mars-size embryos expected in the early asteroid 
belt. They reject the lower disruption criteria as being 
inconsistent with observations (and their initial assumed 
population of 1 km-radius bodies). 

Here, we stress that a population of comparably sized 
bodies will be subject to the lowest possible disruption 
criteria, Q%jj ~ = i- For example, two colliding bodies 
with individual radii of 1 k m a nd density of rock have 
Rci ~ 2 km. Using equation 30 with c* = 5 and p — 0.37 
for small bodies, <%£>,7=i = 5-3 J kg -1 , and the corre- 
sponding V* =1 is 6.5 m s _1 . For a 45-degree impact, the 
value for Q% D rises by almost a factor of 5 (equation 



23) and V* = 14 m s" 1 . Note that this value of Q% D is 



more than a factor of 100 lower than the 5 km s strong 



basalt disruption curve (Figure 12). This example illus- 
trates the need to incorporate self-consistent disruption 
criteria that account for impact velocity and mass ratio 
in order to infer the magnitude of collisional evolution in 
a given population of bodies. 

6.3.3. Application of collision scaling laws in planet 
formation 

To date, all numerical implementations of fragmenta- 
tion during collisional growth of planets assume pure en- 
ergy scaling. That is, p is assumed to be 2/3 and, thus, 
there is no velocity dependence in the catastrophic dis- 
ruption criteria (equation J3]). In contrast, the vast col- 
lection of data in the gravity regime indicate that catas- 
trophic disruption is best fit by nearly pure momentum 
scaling. With nearly linear dependence on the critical ve- 
locity, the criteria for fragmentation m ay vary by orders 
of ma gnitud e during planet formation ( Stewart fc Lein- 
hardt 2009}. Several studies have investigated shifting a 
single reference size-dependent disruption curve by a con- 
stant value that is fixed over the course of the simulation 
to represent stronger or weaker bodies, but none have 
considered a fragmentation criteria that may be variable 
in time and space. 

Furthermore, planet formation models have not in- 
cluded the dependence on the mass ratio of the bodies 
on the disruption criteria. A recent simple treatment of 
the collision parameters required for hit-and-run versus 
merging indicated that planet formation was only slightly 
delaye d compared to simulation s with only merging out- 
comes ( jKokubo fc Genda||2010 ). However, this work did 
not include any treatment of fragmentation. Based on 
our calculation of the region of partial accretion for non- 
grazing impac ts a nd projectile disruption in hit-and-run 
events (Figure 11 ), fragmentation is a critical component 
of the end stage of planet formation. Numerical simula- 
tions cannot assume pure merging or pure hit-and-run. 
The diversity of collision outcomes during the end stage 
of planet formation is described in detail in the compan- 
ion paper, which uses the impact parameters from recent 
TV-body simulations that assumed perfect merging to cal- 
culate the range of collision outco mes predicted by our 
model ( [Stewart fc Leinhardt||2011[ ). 

In this work, the general catastrophic disruption law 
relies upon only two independent mate rial constants {p 
and c*; q g is related through equation 29) and the im- 
pact conditions (mass, mass ratio, impacF angle and ve- 
locity). The material parameters are well constrained. 
The coupling parameter p is tightly constrained by a 



large dataset (Figure |4| to fall close to pure momentum 
scaling (0.33 to 0.37). 

The energy dissipation parameter c* is a measure of 
the physical characteristics of the body. For small bodies 
with varying composition and strength, we find c* = 5 ± 
2. As bodies grow into protoplanets (~ 1000 km), they 
heat internally from the heat of accretion and radioac- 
tive decay. Then, the energy dissipation parameter for 
hydrodynamic planets and planetesimals, c* = 1.9 ±0.3, 
is appropriate. After molten planetesimals recrystallize, 
they will be temporarily stronger until they experience 
sub-catastrophic shattering impact events. Based on 
the strong rock simulations, collisional evolution mod- 
els should test for sensitivity to values of c* up to about 
20. 

Two additional parameters describe the distribution 
of fragments during erosive collisions. /3 is the exponent 
to the size distribution of small fragments and rj is the 
slope of the power law size distribution for the largest 
fragment in the super-catastrophic regime. The values 
for j3 and 77 are also tightly constrained by simulations 
and laboratory experiments, respectively, and may be 
considered, to first order, constant. 

7. CONCLUSIONS 

We present a completely self-consistent set of scaling 
laws to describe the dynamical outcome of a collision be- 
tween any two bodies in the gravity regime. The scaling 
laws rely upon the concept of the interacting mass, the 
fraction of the projectile involved in the collision. Using 
the interacting mass, we derive a general catastrophic 
disruption scaling law. The general forms include two in- 
dependent material parameters: the coupling parameter 
p and the energy dissipation parameter c*. The catas- 
trophic disruption criteria is used to bound the transi- 
tions between collision outcome regimes. The subsequent 
description for the size and velocity of collision fragments 
are derived in closed-form analytic equations with two 
well-constrained parameters. 

With this powerful new tool to describe the dynamical 
outcome of collisions, the physics of collisions in planet 
formation models will have much greater fidelity. Planet 
formation models now have a very small number of free 
parameters to describe collision outcomes (primarily the 
energy dissipation parameter, c*). With a more robust 
physical model for collisions, more detailed factors may 
be studied, such as the evolution of composition dur- 
ing planet formation. Alternatively, other fundamental 
aspects of planet formation may be investigated more 
deeply (e.g., migration) now that the collision model is 
tightly constrained. 

Future work should now bring greater scrutiny to the 
scaling laws used in the strength regime. Although 
this regime has historically been better constrained by 
the abundance of laboratory experiments, self-consistent 
scaling laws also require attention to the dependence of 
fragmentation on the impact velocity, mass ratio, and 
impact angle. 

By fully constraining the dynamics of collisions in the 
gravity regime, this work represents a major advance- 
ment in the robustness of simulations of planet formation 
and the collisional evolution of populations of planets. 
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TABLE A.l 

Summary of primary variables and annotations. 



Symbol 



Definition 



Material parameters 

c* 

I" 

P 

V 

Principal variables 

V,V Z 

Vcsc, V in f 

Q 

Qr 

Q"rd 
M,m 

Rci 
P 
7 
a 

1g 

t 

v 

P 

Geometric variables 

9 

b 

frcrit 

R 
r 
D 

l/(2r) 

Superscripts 



Subscripts 

targ 

p 

tot 

interact 
a 

-1=1 

lr 

sir 

rem 

core 

Constants 



Pi 
G 



Head-on equal-mass disruption energy in units of specific gravitational binding energy 
Velocity exponent in coupling parameter 
Slope of fragment size distribution 

Exponent in power-law fragment distribution in super-catastrophic regime 
Impact velocity 

Escape velocity, velocity at infinity 
Specific impact energy 

Specific impact energy for the collision in center of mass frame 

Catastrophic disruption criteria - specific impact energy to disperse half the total mass 
Mass 

Radius of total mass in a body with density 1000 kg m -3 
Reduced mass 

Projectile-to-target mass ratio 

Mass fraction of the projectile that intersects the target 
Coefficient of gravity term in general formula for Q* RD 
Accretion efficiency 
Velocity of fragments 
Density 

Impact angle (0 is head-on) 
Impact parameter sin(f?) 

Critical impact parameter denoting transition from non-grazing to grazing collision 
Target radius 
Projectile radius 
Diameter 

Fraction of projectile diameter that overlaps with target 

Value at the catastrophic disruption threshold 
Oblique impact 

Reverse impact onto the projectile in the hit-and-run regime 

Target 

Projectile 

Target + projectile 

Interacting fraction 

Interacting projectile fraction 

Equal-mass collision 

Largest remnant 

Second largest remnant 

Fragments smaller than the largest remnant 
Core fraction of a differentiated body 

Density of 1000 kg m -3 
Gravitational constant 



A description of variables and annot atio ns is given in Table A.l The compilation of simulation data on catastrophic 
disruption presented in Figures [4] and 12 is s umm arized in Table A. 2 A summary of all the PKDGRAV simulations 
conducted in this study is presented in Table |A.3[ 

8.1. A general formulation for collision outcomes 

We summarize the sequence of logic that should be applied to determine the dynamical outcome of any collision in 
the gravity regime using our analytic model. First, we identify the boundaries of the major collision regimes: 

1. For a given collision scenario (M p , M targ , 6, Vi, and R p and i?tar g from the bulk densities of the bodies), calculate 
the interacting mass fraction of the projectile, TOj n t er act = o:M p (equation 11). 



If Vi < V^ sc (equation 53 ) , then the impact is in the perfect merging regime. 

Calculate the critical impact parameter 6 cr it for the collision (equation [6]) . If b < b CT it, then it is a non-grazing 
impact, else it is a grazing impact. 



Calculate the catastrophic disruption criteria, Q'^jj, and corresponding critical impact velocity, V 
specific impact scenario: 



for the 
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TABLE A. 2 

Sources and description of catastrophic disruption data presented in Figures [4] and |T2j Filled and line symbols indicate 

HEAD-ON COLLISIONS (6 = 0). OPEN SYMBOLS INDICATE OBLIQUE IMPACTS: b = 0.15,0.3 FOR OPEN STAR, 0.35 < b < 0.9 FOR A, b = 0.5 FOR 

□ AND 0, AND b = 0.71 FOR XI, HOURGLASS, HEXAGON, ®, AND E3. 



Symbol 



Target Description 



Source 



Weak targets 

A, A 

+ 

if, open star 



hexagon 
Strong targets 
T 

X 

hourglass 

Hydrodynamic planets 



x 



10-km PKDGRAV rubble piles 
1 to 50-km PKDGRAV rubble piles 
1-km PKDGRAV rubble piles 
weak 2 to 50-km basalt 
hydrodynamic 2 to 50-km basalt 
50-km ice 



1 to 100-km radius basalt 

2 to 50-km basalt 

0.3 to 100-km microporous rock (pumice) 
0.3 to 100-km basalt 
10-km macroporous basalt 
10-km basalt 



2.2 Mercury-mass bodies, differentiated rock and iron 
0.4 and 7 Earth-mass bodies, pure rock 
1 to 10 Earth-mass bodies, differentiated rock and iron 
0.5 to 5 Earth-mass bodies, differentiated water and rock 



This work 

Stewart fe Leinhardt| |2009| 
,einhardt et al.|j20flo) 
■cinhardt & Stewart (2009) 
.cinhardt & Stewart (2009) 



Leinhardt fc Stewart| ( |2009| 



Benz fc Asphaug|l|l999| 



Leinhardt & Stewart 
Jutzi et aL| | |2010| 
Jutzi et aiqj20l"0" 
Benz| ||2000| 



1 2009 I 



Benz (2000) 



|Benz et al.| (|2007[) 
Marcus et al.| {2009 
Marcus etal"] |20(39 i 
Marcus et^DpOlO 



(a) Calculate Rci from the total mass and density of 1000 kg m 3 . 

(b) Calculate the principal disruption value for an equivalent equal-mass collision at Rex, Q*rd.^=i (equation 
28), and its corresponding critical impact velocity, V* =1 (equation 30), using the material parameter c*. 



(c) Calculate the reduced mass, pL, and the reduced mass using the interacting mass, fi a (equation 12). 

(d) Calculate the disruption criteria, Q*f> n , and critical impact velocity V* , for a head-on impact with the 
desired mass ratio, 7 using equations |23| and |22| and the material parameter fl. 

(e) The value for the disruption energy, Qr D , and critical impact velocity V R * D , for the desired impact angle 
are found using equations |15| and |16| 

Calculate the value for Qr required for onset of erosion, M\ T — -Mtarg, using the value of Qr D and the universal 
law for the mass of the largest remnant (equation [5]). From this Qr, derive the impact velocity for the onset of 
target erosion, Ve ros i on , from equation [T] 



6. For grazing impacts (b > & cr it), the hit-and-run regime is bounded by V" ' sc < Vi < V c 



erosion • 



Note that the 



graze- a nd-merge regime is a subset of this range, but it is not explicitly defined in this work (see Genda et al. 
2011b|. 



7. Calculate the value for Qr required for the onset of super-catastrophic disruption, M\c = 0.1Af tot , using the 
value of Q'rd and the universal law for the mass of the largest remnant (equation [5]) . From this Qr, derive the 
impact velocity for the onset of super-catastrophic disruption, V^ uperca t, from equation [TJ 

8. For all impact angles, the target is eroded when Vt > V^rosion and the impact is in the erosion regime. 

9. For all impact angles, the impact is in the super- catastrophic disruption regime when Vi > Kj U p Crcat . 

10. For non-grazing events and V^ sc < Vi < Ks Uper catj the impact is in the disruption regime and the universal law 
for the mass of the largest remnant applies. The partial accretion regime is bounded by V^' sc < Vi < Vision- 



11. For grazing events and V c 

the mass of the largest remnant applies only for Mi r < M tar g 



< v<<. m < Vi < tgupercat) th° impact is in the disruption regime and the universal law for 



12. In the hit-and-run regime, calculate the critical disruption energy for the reverse impact onto the projectile, Q*r D 



as described in £4.2 



and its corresponding * . Use the equation [5] or [44] to determine the largest remnant after 



M p . 



disruption of the total mass involved in the reverse impact, Mi n t era ct 

In the merging regime, mass and momentum are conserved. 
In the disruption regime: 

1. Determine the mass of the largest remnant M\ x from the universal law (equation [5]) using Qr and Qrjj- 
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2. Determine the mass of the second largest remnant M s i r using equation [37| with ft = 2.85, N\ = 1, and N2 — 2. 



The size distribution of the tail of smaller fragments is described by equation 31 



3. For b = 0, assume that the largest remnant obtains the velocity of the center of mass; for b > 0.7 assume that 
the largest remnant maintains Vtarg- For < b < 0.7, the largest remnant velocity is some quasi-linear function 
of b. 

4. The mass-velocity distribution of the smaller fragments with respect to the largest remnant is given by equation 



In the super- catastrophic disruption regime: 



1. De termine the mass of the largest remnant M\ t from the power law (equation 44) using Qji and Q% D (equation 



15). 



2. The size and velocity distribution of the fragments with respect to the largest remnant are the same as in the 
disruption regime. 

In the hit-and-run regime: 

1. The mass of the largest remnant M\ r is approximately equal to the target mass M targ . 

2. The mass of the second largest remnant M srr is estimated using the universal law and disruption criteria for the 
reverse impact on the projectile, Q^* D - 

3. When the projectile is disrupted, the size and velocity distribution of the fragments are described as in the 
disruption regime with respect to the largest remnant from the projectile. 

4. In the special case of 7 ~ 1, the onset of erosion occurs simultaneously in both bodies and M\ r ~ M s \ r . Then, 
use Ni = 2 and N2 = 4 to calculate the size distribution. One can assume that the fragments from both the 
projectile and target have identical size and velocity distributions with respect to their body of origin. 

In the disruption and super-catastrophic disruption regimes, a differentiated target may change its bulk composition 
by stripping off a portion of the mantle. Bulk compositional changes may be tracked using the results from Marcus| 
et al. (2010b I. They found that the core mass fraction after a disruptive collision falls between two idealized models, 
and we suggest using an average of these two results: 



1. Model 1 - Cores always merge: Given the original M corc ,targ and M< 



min(Mi r ,M cor0itarg + M corc , p ). 



core, p ; 



the post-impact core is M c , 



2. Model 2 - Cores only merge on accretion: When M\ r > M tar g, M coro = M corc ,targ + min(Af coreiP , M\ r — M targ ). 
When Mi r < M targ , assume that none of the projectile accretes and the mantle is stripped first. Then, M corc = 
min(M corc , ta rg,Mi r ). 

In hit-and-run events with projectile disruption, the same relations may be used to estimate the bulk changes in 
compositi on for the projectile. 



Finally, Housen & Holsapple (2011 1 provide scaling laws for the gravitationally escaping ejecta from the target in the 
impact cratering regime. The impact cratering regime occurs at the low velocity en d of the disruption regime, w hen 
Vi > V esc , Al p « Mt ar g, and M lr ~ M targ . Based on many laboratory experiments, Housen & Holsapple (20111) find 
that approxim ately 0.01Af p achieves escape velocity in cratering events of Vi ~ V esc (see their Figure lb). In addition, 
Svetov (20111 provides empirical equations for ejected material in cratering collisions on self-gravitating bodies. 



Collision Scaling Laws 



TABLE A. 3 
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